H 7 2 


28 8 1 9 


« 


Mobil 




pX, 


FINAL REPORT 

NEUTRON DIE-AWAY EXPERIMENT FOR REMOTE ANALYSIS 
OF THE SURFACE OF THE MOON AND THE PLANETS - 
PHASE I I I 


MOBIL RESEARCH AND DEVELOPMENT CORPORATION 


FIELD RESEARCH LABORATORY 
DALLAS, TEXAS 


FINAL REPORT 

NEUTRON DIE-AWAY EXPERIMENT FOR REMOTE ANALYSIS OF 
THE SURFACE OF THE MOON AND THE PLANETS - PHASE I I I 


Period from 1 June 1970 to 30 November 1971 


Prepared by: W. R. Mills and L. S. Allen 

Mobil Research £• Development Corporation 

Field Research Laboratory 

3600 Duncanville Road 

Dallas, Texas 75211 

214-331-6531 


For: Robert P. Bryson 

Apollo Lunar Exploration Office 
Code MAL 

National Aeronautics and Space Administration 
Washington, D. C. 20546 




TABLE OF CONTENTS 


Paqe 


FOREWORD 


ABSTRACT 


LIST OF TABLES 


LIST OF FIGURES 


INTRODUCTION 


II. DATA ANALYSIS AND PREDICTION OF VARIANCE 

A. Least Squares Analysis 

B. Prediction Analysis 

C. Systematic Errors 

III. EPITHERMAL DIE-AWAY 

A. Functional Relationships Among Variables 

B. Experimental Calibration Studies 

C. Theoretical Calculations Related to Experiments 

D. Status of Analysis of Epithermal Die-Away Data 

IV. THERMAL DIE-AWAY 

A. Background Discussion 

B. Determination of S Values for Rock Models 

a 

C. New Thermal Neutron Die-Away Measurements on 
Large Model s 

D. Method of Data Analysis 

E. New Capture Gamma-Ray Die-Away Measurements 


V. ANALYSIS OF VENUSIAN ATMOSPHERE AND SURFACE BY 
LIGHTWEIGHT CPNE PROBE 



Page 


VI. MISCELLANEOUS TOPICS 99 

A. Neutron Output Monitor 99 

■? 

B. Flight-Qualified He Detector 103 

VII, SUMMARY 105 

LIST OF SYMBOLS IO7 

REFERENCES IO9 

APPENDICES 

TABLES 

FIGURES 


-i i- 



FOREWORD 
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Co-investigator 
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ABSTRACT 


Continuing work on the two die-away measurements proposed to be 
made in the Combined Pulsed Neutron Experiment (CPNE) for analysis of 
lunar and planetary surfaces is described. The present report documents 
research done during Phase III. 

A general exposition of data analysis by the least-squares method 
and the related problem of the prediction of variance is given. This 
exposition is based on a book. Prediction Analysis, by J. R. Wolberg. 
Prediction analysis is a technique that is useful for estimating ex- 
pected variances of measured variables before the data are available. 

Both analysis of variance and prediction of variance are important in 
all of the measurements carried out in the CPNE, and are applied 
specifically in this report to the epithermal and thermal neutron die- 
away experiments. 

A data analysis procedure for epithermal die-away data has been 
formulated. In order to facilitate the analysis, the number of inde- 
pendent material variables has been reduced to two: the hydrogen density 

and an effective oxygen density, the latter being determined uniquely 
from the non-hydrogenous elemental composition. Justification for this 
reduction in the number of variables is based on a set of 27 new theo- 
retical calculations. Work is described related to experimental calibra- 
tion of the epithermal die-away measurement. Certain difficulties were 
encountered and necessitated an extensive theoretical study of the new 
experimental results. These studies lead us to believe that ultimately 



the best calibration procedure will involve a combination of theory and 
experiment. An interim data analysis technique based solely on theo- 
retical calculations seems to be adequate and will be used for future 
CPNE field tests. 

A method of analyzing thermal neutron and capture gamma-ray die- 
away data has been devised. In the process of investigating this method 
it was found that two basic difficulties existed with previous thermal 
die-away measurements. The measurements on the large rock models were 
distorted by late returning neutrons, and the values of the macroscopic 
thermal neutron absorption cross section deduced from chemical analyses 
were in error. Both of these difficulties have been corrected by new 
measurements which are described. The new method of analysis permits a 
complete empirical calibration to be carried out, and this method will 
be used on the thermal die-away data obtained in future CPNE field tests 

A brief study has been made on the expected performance of the 
epithermal die-away measurement in analyzing the atmosphere of Venus 
for hydrogen. Although no numerical results have been obtained, it is 
proposed that the study could form the basis of more extensive investi- 
gations in the future. 

A new method of remotely monitoring the neutron output of the 

generator used in the CPNE is proposed and investigated. This method 

3 3 

is based on the detection of He recoils in the He/Cd detector used 
for the epithermal die-away measurement, and would require the addition 
of only a pulse discriminator and scaler to be implemented. 
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I . INTRODUCTION 


During Phase I of work carried out under Contract NASW-l435^ 

"Neutron Die-Away Experiment for Lunar and Planetary Surface Analysis", 
the feasibility of performing thermal and epithermal neutron die-away 
measurements on remote planetary surfaces was i nvesti gated , ^ ^ ^ It was 

demonstrated that such measurements are feasible under reasonable ex- 
perimental conditions. During Phase II of Contract NASW-1435, the work 
was extended to include some problems of data analysis, possible con- 
flicts with other parts of the Combined Pulsed Neutron Experiment (CPNE), 

and a demonstration of the feasibility of the CPNE concept using a single 

( 2 ) 

experimental probe and data acquisition system.' ' 

On 19 and 20 November 1969^ a meeting of personnel associated with 

the CPNE was held at Goddard Space Flight Center to discuss the status 

of the work and future plans. A list of "action items" evolved from 

that meeting, and this list formed the broad basis for proposed work 

to be done during the next phase. The complete list of items has been 

(3) 

documented by Reed and handler. ' The most significant and critical 
points are related to the problem of choosing a final probe configuration 
for the experiment. The configuration must be selected as a result of 
a reasonable compromise among several parameters, including weight, 
volume, power, experiment running time, and quality of the results ob- 
tained. The configuration is also influenced by other factors, such as 
the radiation field from prospective RTG electrical power sources, the 
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ambient cosmi c-ray- induced neutron and gamma-ray fields present at the 
surface of a planet, and any unusual requirements which the CPNE would 
impose on the host vehicle. 

The question of probe weight became the prime consideration as a 
result of a presentation made by R. L. Caldwell and W. R. Mills of 
Mobil and J. H. Reed of IITRI to planetology personnel at NASA Head- 
quarters on 22 October 1970. From this presentation and associated 
discussions, it became clear that the CPNE may occupy a unique position 
in considerations regarding surface missions to Venus. At present an 
experiment utilizing neutron techniques seems to be the only one which 
is feasible for elemental analysis. This is due to the very high sur- 
face pressure (100-200 atm) believed to exist on Venus. However, due 
to the small science payload weight which will likely be allotted to 
geochemical experiments, it is absolutely essential to reduce the weight 
of the CPNE to the minimum possible which still allows useful measurements 
to be carried out. 

In its simplest and most relevant terms, the question of the 

probe configuration and weight is as follows: If the values of n physical 

variables, p^, ... p^, are either directly measured or deduced from the 

2 2 

various parts of the CPNE, what are the expected variances, a , ... ct , 

Pi Pn 

of the values for a given configuration, weight, and environmental con- 
ditions? For a given probe configuration and weight, neutron output, and 
data collection time, the standard deviation a for the i ' th variable 

Pi 

is one measure of the sensitivity (or threshold for measurement) of the 
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CPNE for that variable. A more realistic estimate of the true sensi- 
tivity is three or four times a , but the standard deviation (or the 
variance) is a convenient measure of sensitivity and this convention was 
agreed upon at the GSFC meeting. The variances of the variables are 
also convenient indicators of the influence of the RTG radiation field 

and/or the ambient neutron and gamma-ray fields produced by cosmic rays. 

2 

These effects can only result in increasing the a , thereby increasing 

^i 

the uncertainty on a given measurement and increasing the measurement 
threshold for the variables. 

From the above paragraph, it is evident that realistic and defensible 
predicted estimates of the variances of all measured variables in the 
CPNE are essential in order to objectively decide on the most desirable 
ultimate probe configuration for any given mission. For this reason, a 
large portion of our effort during Phase III has been addressed to the 
problem of determining the variances of the variables which are measured 
in the die-away portions of the CPNE. In order to accomplish this, the 
procedures for analyzing epithermal and thermal neutron die-away data 
must be specifically worked out and the relation between the data analysis 
and calibration procedures carefully studied. The bulk of the present 
report deals with this problem. 
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II. DATA ANALYSIS AND PREDICTION OF VARIANCE 


A. Least Squares Analysis 

The method which we will use to obtain or predict the variances 
of variables deduced from die-away measurements is that given by 

( 4 ) 

Wolberg. One assumes that the data will be analyzed by the least 
squares technique. The solution of the most general case of least 
squares analysis is given by Wolberg, including least squares estimates 
of variances and covariances of adjustable parameters. Predicted values 
of the parameter variances and covariances are obtained by using certain 
assumptions and the formal results of the least squares analysis. This 
procedure will be sketched briefly below in order to provide a self- 
contained basis for the application of the method to the die-away data. 
Further details can be found in the very readable book by Wolberg. 

Symbols used in the discussion below are given in the list at the end 
of this report. 

Suppose that we obtain experimental values of a dependent variable 
Y. and m independent variables . (i=1,n; J=1,m) for n observations or 
data points. We assume that the true values of these variables are Tlj 
and I j . . If systematic errors are absent and if the experiment consisting 
of the n observations were repeated a large number of times, in the limit 
the mean value of Y. approaches 1]. and the mean value of X ^ . approaches 


5jr 


As a first step in least squares analysis, a functional relationship 
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between the variables is assumed: 




(1) 


The function f is chosen on the basis of whatever a priori informa- 
tion is avail able^ and contains the parameters (k=1,p). The oTj^'s 
may or may not be directly related to physical variables. The purpose 
of least squares analysis is to use the available data to find "best 
estimates" of the q'j^. These are the least squares values aj^. If 
systematic errors are absent, and if the experiment is repeated many 
times, the average values of the aj^ will approach the For a given 

set of experimental data, calculated or least-squares adjusted values 
of the dependent and independent variables are assumed to be related 
by the functional form in Eq. (l) with qCi^ = aj^; k=1,p: 


y. = f (x, . . . X . ; 
^1 ' 1 r mi 


; a 


r • 




( 2 ) 


In Eq. (2), y. and x.. are the adjusted values of the variables. 

(4) 

Figure 1, taken from Wolberg shows the relation between the 
observed, adjusted, and true values of the variables for the case of 
one independent variable. The residuals and errors are defined by 


R = Y. - y. 
y. I ' I 


(3) 
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R = X. . - X, , 
X. . J I J I 

j I ■> 




E = Y. - 11. 

y. I 


( 5 ) 


E = X.. - 

X . . J I II 

ji J 


( 6 ) 


In general, E and E can only be estimated, 
y • X • • 

' J ' 

The weighted sum of squares of residuals is 


m 


= Z(w R^+ Zw R^) 


1=1 


. , X. . X. . 

J=1 JI JI 


(7) 


where w and w are suitably chosen weighting factors. The 

y • X • • 

I J ' 

method of least squares consists of the determination of the values 
a. which minimize S. In many problems the residuals R can be 

K X • • 

J I 

neglected. In the general case they must be retained in S. The 

form in Eq. (7) is correct provided that the errors E are uncor- 

^i 

related and the errors E are uncorrelated. 

X. . 

J I 

Wolberg gives a solution of the minimization problem in an itera- 
tive form. He defines quantities where 


\ ' =ko ■ =k ••• 


(8) 


and aj^^ is an initial guess for the value of aj^. Letting A be a 1 x p 
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row vector, the solution is given as 


A = 



V, 


where C ^ is the inverse of a p x p matrix and V is a p x 
vector. The elements of C and V are given by 





) 


where the quantities F 


i 

a 


k 



and L. 

I 


are defined by 



9f (x, . .X . ; a,, . . .a ) 
1 r mi V o' 

da, 


X . . “X ..,..«x • 
1 1 1 r mi 


F 


i 

0 


Y. 

I 



.X . 
mi 


10 ’ 



(9) 

1 column 


( 10 ) 


( 11 ) 


m I 


(13) 
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/c* \2 

m (F ) 

tn ^ • 


L. - — . I 

y. j=l X.. 

U J j I 


1i mi V p 


X, .=X . .X .=X . (15) 

1 1 1 r mi mi 


Initial guesses are made for aj^^, the are calculated from Eq. 

(9)j and new guesses a^^^ - Aj^ are thus obtained for the next iteration. 
The process is continued until a convergence criterion is met. With 
some non-linear problems, it may happen that a reasonable convergence 
criterion cannot be met. Wolberg discusses the steps one then takes 
to try to obtain convergence. 

The second part of the least squares analysis. is the estimate 
of the variances and covariances for the parameters. The variance 
of aj^ is defined by 


‘’a, = Gk ' 0 ’ 


where the bar signifies the mean value of the squared quantity which 

is obtained when the experiment is repeated a large number of times. 

The covariance between a, and a, is 

k 1 


^kl ” ^kl 


''a, "a. = [ • "k> <=1 ■ "l* ] ’ 


'k °1 
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where p, , is the correlation coefficient between a, and a,. Under 
kl k 1 

the assumptions that the errors are not large and that .the quantities 
E. are uncorrelated, where 

I ’ 


m 


E. = E + 

I y j X . X . . 

’'i J=1 J Ji 


J) F* E (i=1,...,n) 


Wolberg shows that 


2 S - -1 

*^a. n - p kk 


(18) 


and 


o, , « 

kl n - p 


Cj', k ^ 1 


(19) 


where S is the weighted sum of the residuals for the converged values 

of the Equations (18) and (19) are approximate because the 

data of a single experiment can yield only estimates of the variances 

2 

and covariances. Unbiased estimates of cr and are denoted by 
2 

s and s , respectively, and are given by 
®k 


2 

s = 


-1 


n - p 


'kk 


( 20 ) 


’kl 


n - p 


'kl 


( 21 ) 
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In the applications with which we are concerned (i.e, die-away 

2 

measurements), we are not interested in the values of a and ct, , 

.3 1 K I 

2 ^ 

as such. Instead, we are interested in o when y is determined from 

y 

a measurement on an unknown material, assuming that the have 

been determined in a previous calibration. Thus, we have 


y = f(x^,... x^; a^,... a^), (22) 


2 

and the measured variables x. have variances a (and possibly non- 

■* 

zero covariances a ). The variance of y is determined as follows. 

X.X.' ’ 

' J 

From Eq. (22), 


m p 



(23) 


By the definition of variance. 


a 


2 

y 



(24) 


Squaring Eq, (23) and averaging, we have 
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m p 



(25) 


If we assume no correlation between the Ax. and Aaj^, which is a 
reasonable assumption, then 


Z f (t)(i,) ^ t t (i,)(i,) 


i=l j=1 ' J J k=l 1=1 


m m 


P P 


I I °x.x. * Z Z %, °kl 


i=1 j=1 ' J ' J k=1 1=1 


'k °1 


2 2 

In Eq. (27), a = a , o. , = <7 , and the partial derivatives F 

^ x.x. X. ^ kk a, X. 

I I I k I 

and F are to be evaluated for the measured values of x,,..., x 
a, ^ ' m 

^ 2 
on the unknown material. An estimate of cr is given by 


mm P P 

X. - Z Z 

i=l J=1 ' J ' J k=1 1=1 ' 


= Z Z f f ^ 

Li Lt X. X. X. 


where s i s an estimate of the covariance between the independent 
x.x. 

I J 

variable measurements, and s., is given by Eq. (21). Also, s = 

2.2 ' ' 
s and s, , = s 
X. kk a, 

I k 
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B. Prediction Analysis 

We now wish to answer the following question: Assuming the data 

are subjected to least squares analysis, can we predict the values of 
2 2 

a , CT. ,, and o in Eqs. (18) (19)^ and (27), respectively? These 

^ ^2 2 
predicted values are denoted by o , and cr . Similarly, the 


predicted values of any other quantities will be shown with a caret 
(4) 

on top. Wolberg adopts the definition that the predicted value 
of any quantity is the mean value of the unbiased estimate of that 
quantity which would be obtained if the experiment were repeated many 
times. Thus, 



The quantity S has the "chi -squared" distribution, and therefore has a 

A 

mean value equal to the number of degrees of freedom. Thus, S = n-p, 
and we have 



( 32 ) 
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= Q 


2 

a = 

y 


tn m 


2 z 

i=i j=i 


A /\ /\ 

F s + 

X. X.x. 

I J » J 


z z 

k=1 1=1 


A A 

F F 


A \-1 


(C,) 


(33) 


(34) 


Equations (32) - (34) are the basic equations of prediction analysis. 

In order to carry out prediction analysis calculations, one must 
solve two main problems: 

(i) Obtain a suitable functional form f(x^,... x^; 0 ^^) 

from which partial derivatives can be evaluated, and 
(ii) calculate the predicted elements of the matrix C, from 
which the predicted inverse (^) ^ is obtained. 

The predicted elements of C are given from Eq. (10) by 



z 

• 




L. 


(35) 


C. Systematic Errors 

The treatment of least squares analysis, and the subsequent pre- 
diction analysis, was based on the assumption that only random errors 
are present. If the experiment is repeated a large number of times, 
the average values of all quantities would approach their true values 


Y. 

I 
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a, - O', 

k k 

If systematic errors are present, some or all of the above averages 

(4) 

will not approach the true values. Wolberg' classifies two types 
of systematic error: 

Type 1 : Errors which cause the measured values of a variable 

to not approach the true value. 

Type 2 : Errors which cause the average aj^ to differ from the 

even though type 7 errors are absent. 

Type 1 errors would normally be expected to arise from purely experi- 
mental sources. Type 2 errors would probably arise most often in the 
data analysis process. Examples of type 2 errors are, 1) using an in- 
appropriate function in Eq. (1), and 2) failing to include one or more 
independent variables in the analysis. 

We will not go further into the general problem of systematic errors 
here, in effect, one must identify all the possible sources of systematic 
errors, determine how they enter the analysis, and estimate the magnitude 


of the errors. 
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III. EPITHERMAL DIE-AWAY 

A. Functional Relationships Among Variables 

At the conclusion of the work done on epithermal die-away during 
Phase II, two principal problems remained. These were (1) the deriva- 
tion of a convenient and workable relationship among experimental data 
obtained in an epithermal die-away measurement and the material neutron 
constants which govern the die-away, and (2) the development of a sat- 
isfactory calibration procedure. This section is devoted to the first 
problem. 

We first need some measure of an observed epithermal die-away curve. 
Experimental data for the four large models at Mobil were given in Fig. 

17 of Reference 2. For ready reference, these data are reproduced here 
In Fig. 2. As a measure of the die-away rate, we have chosen to use 
the approximate negative logarithmic slope at a delay time of 50 |j,sec. 

If R(t) is the actual detector response as a function of delay time t, 
then the negative logarithmic slope at 50 |j,sec is 




t=50 


For actual data, we will approximate InR(t) by a straight line over 
the range 40 to 60 psec; 



40 — t — 60 psec. 


InR(t) 


-yt + b. 
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The letters l.s. over the equals sign indicate that the parameters y 
and b will be determined by least squares fit; thus 


y Ri 




( 36 ) 


The choice of the negative logarithmic slope, as determined over the 
range 40 to 60 psec, for the experimental parameter does not, of 
course, make use of all the epithermal die-away data obtained in a 
measurement. Use of the complete set of data points would lead to a 
very difficult non-linear least squares problem in order to deduce 
estimates of material neutron constants from the data. The development 
in this section leads to a linear least squares analysis and this must 
suffice for the present. 

We now need to decide what the appropriate variables are which 

describe the physical and chemical properties of the material we are 

trying to analyze. To do this, consider the material parameters which 

enter the Boltzmann transport equation for epithermal neutrons. We 

previously used the Greul i ng-Goertzel scattering kernel approximation 

and the diffusion approximation to obtain an equation for the neutron 

flux 0(r,z,u, t), where r and z are cylindrical coordinates, u is neutron 

( 2 ) 

lethargy, and t is time. The equation, good only for neutrons with 
energy greater than about 1 ev and for constant scattering cross section, 
is 


— e 


2 ^ 
dt 


0 - D 


V^0 


IT “ 

5 J e 0 du' + S(r,z,u,t) 

Y 0 


( 37 ) 
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2 

In Eq. {yj) , Vq is the speed of source neutrons, V is the Laplacian 
operator in cylindrical coordinates, and S is the source term. The 
physical parameters S, D, §, and y a*"® given below for a single element 


E = 


0.602 

M 


P CT 


( 38 ) 


D = 


3(1- Im) S 


( 39 ) 


^ 2M \M-1/ 


(40) 




In these equations, 

M = mass of scattering nucleus relative to the neutron mass, 
p = bulk density (gm/cc), 

a = microscopic scattering cross section (barns). 

For mixtures of elements. 


S 


2 


0.602 


w. 

I 


CT. 

I 


m 
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1 y 0.602 

Z L M. 

. I 


p w. a. §. 

^ I I ^1 


m 


Y = 


1 

S 5 


I 


0.602 

M. 


P «5 a. Y, 


(45) 


where w. Is the weight fraction of the i ' th type element, and the 
other subscripted quantities refer also to the i ' th element. For 
a given suite of elements present, all of the quantities in Eqs. 

(42) - (45) are known nuclear constants except p and the w. 's. 

These latter quantities are the proper material variables. 

Since we hope and expect that the epithermal die-away measure- 
ment will be sensitive primarily to hydrogen, it is convenient to 
separate its contribution in Eq. (42) - (45). If we define 


r. “ 

I 


0.602 

M. 


ai , 


(46) 


the epithermal neutron parameters may be written as follows; 
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M2.64 w + 
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(47) 
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- p[;i2.64„„+3 Z r. w, (I - , 

I ^ 


m 


§ = 


12.64 w + y r. w. 

H ^4 I I I 


12.64 + r r. w. 


(49) 


Y = 


12.64 w,, + ) r. w. §. V. 

H 4j I I I 'i 

I 

/ 

12.64 r. w. 


(50) 


The primes on the sums mean that H is omitted. The form of Eqs. 
(4?)- ( 50 ) suggests that a simplification can be made by using the 
elemental densities: 


P, = P w. . 


With this change, we have 

t 

E = 12.64 P„ + ^ r. p. 

i 


/ 

D-' = 12.64 3 2 '"r P| (' - 3^), 


(51) 


(52) 


(53) 
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12.64 r. p. C. 

5 - '■ , (54) 

/ 

12.64 P„ +X ‘‘i P. 
i 

/ 

12,64 P„ ^ ^ r, p, ?, V, 

Y - (55) 

12.64 p^ +[[ r. p, 5, 

i 

and the simplest complete set of material variables consists of the 
elemental densities p.. 

As stated above, we expect the epithermal die-away measurement 

to be sensitive primarily to hydrogen. This suggests that we try 

further simplifying the problem by combining the effects of all the 

other elements present into a single equivalent density. If we call 

this variable p then we have a problem in only two material 

non-n 

variables, p,, and p The question is whether this reduction 

H non-n 

can be done with sufficient accuracy and how is p to be calcu- 

' ^non-H 

lated from a given set of p.. If the number of material variables 
can be reduced to two, it will clearly simplify the data analysis 
problem considerably. 

To look at ways of obtaining p ,,, it is convenient to re- 
' ^ '^non-H 

write Eq. (53)* A good approximation for M » 1 is 
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§ 



Since the primed sum in Eq, (53) does not include hydrogen, and 
in practical cases would include only values of M > 12, we drop 
the 2/3 in the denominator. Thus, 


1 I 

3M ~ 3 ^ 


2 1 

For oxygen, = O.Okiy and J § ~ 0.0400. Thus we have a maximum 

2 

error of about 0.2% in the factor ^ , and this is in the 

3M 

parameter 0 which probably does not play a major role in the epi- 
thermal die-away measurement. With the above approximation, we 
can wri te 


/ / 

d“^ » 12.64 Ph + 3 Yi ’’i Pi ~Y. '*i Pi ^i ' 

i i 

From a given set of values of p. (excluding hydrogen) there 
are literally an infinite number of ways in which one could cal- 
culate a p A method based on preserving the correct value of 

at least one material neutron parameter seems the best to try. It 
will probably be impossible to preserve the correct value for all 
neutron parameters. Note that in Eqs. (52), (54), (55), and (56), 
three different primed sums appear; 
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I ''i Pi ' I '■i Pi ^i ’ p"‘‘ Z Pi Pi ^i ''i 


Our approach to obtaining a p ,, is to define effective ele- 

^ '^non-H 

mental densities, relative to a reference element, which will make 
one of the above primed sums correct. Arbitrarily we choose 
oxygen as the reference element and define the following: 


'‘o Pe 


^ Z' Pi Pi 


( 57 ) 


Pfl ?0 P? 5 Pi Pi 5| 


( 58 ) 


Pfl ?0 Yo Py = Z Pi Pi ^i 


( 59 ) 


Thus, 


= r Pi 


( 60 ) 


Y' *^1 ^i 

Li r„ 


Pt 


0 =0 


(61) 
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p 


y 



^ Vj 
•■o ^0 ^0 


( 62 ) 


Use of one of the above three "effective oxygen densities" will 
make the value of T,, or §Sy correct, but not more than one of 
the three will be correct. The problem now is to decide which, if 
any, of the three effective oxygen densities can be used in epi- 
thermal die-away data analysis. Table 1 gives values of the im- 
portant neutron-nucleus constants. The values of a. are taken from 

( 6 ) 

the compilation by Hughes and Schwartz. 

We now have selected variables y, p.., and p ,, (where p ,, 

H non-H non-H 

is either p^, p^, or p^), and must find a suitable functional form 
relating them: 


" '(‘’h’ Pnon-H’ V 

We should try the three expressions for in Eq. (63) and 

select the one which gives best agreement with some "ground truth" 

information. We would still have to decide whether the best way of 

calculating P^^^^ ^ is good enough for our purposes. Since complete 

experimental epithermal die-away calibration data are not available, 

we must use the results of theoretical calculations as our ground 

truth. Also, theoretical calculations provide a much wider range 

of the variables over which to test the choices for p ^ and the 

non-H 
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functional form selected in Eq. (63). To this end we have calculated 
TJ new theoretical epithermal die-away curves. These calcuiations 
are based on the use of Eqs (6.3), (6.19), and (6.21) in Reference 2. 
Each combination possible among the following was calculated: 


Elemental composition: basalt, granite, dunite 

p = 1.5, 2, 2.5 gm/cc 
w„ = 0, 0.0015, 0.003 

n 


Elemental compositions used for the three rock types are those given in 
Table 14 of Reference 2, Table 2 of the present report gives a summary 
of the pertinent parameters. Values of y(in psec ') were calculated 
from the formula ^ In , where R(40) and R(60) are the relative 

theoretical epithermal responses at 40 and 60 psec, respectively. 

The simplest analytic form for Eq. (63) is a linear one: 


y = 


®2 ®3 ^non-H 


(64) 


Sample plots of the results given in Table 2 indicate that such a linear 
relation is inadequate. We are thus led to try a general quadratic 
form: 


y = 


^ =2 Ph 


^ ^3 ‘'non-H ^4 "h 


Pu P u + 

5 H non-H 6 


non-H 


(65) 
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Considering the variables y, p,,. and p ,, as mathematical entities, 
3 u ^non-H ’ 

we could write, instead of Eqs. ( 63 ) and ( 65 ), 


^non-H’ 


( 66 ) 


and 


Pu = + a-y + a.p 


lP u + 3 y,y + ^rY P u + ^£iP u 
3 non-H 4 5 non-H 6 non-H 


( 67 ) 


The advantage of these forms is that p^ (and its variance) is the 

n 

quantity we want to determine for an unknown, and p is explicitly 

given by Eq. (67). It is thus worthwhile to try both Eqs. (65) and 

(67) on the results given in Table 2 . 

The expressions in Eqs. (65)' and (67) were used in a least squares 

fit to the theoretical data in Table 2 . For p ,,, the values p^, 

non-H 2. 

p^, and p^ were used in turn. The sum of squares of residuals was 
taken as a measure of the goodness of fit, although this sum has no 
statistical significance here. For both the forms in Eqs. (65) and 
(67), the use of p^ led to the smallest sum of squares of residuals. 

In order to compare Eqs. (65) and (67), we solved Eq. (65) for p and 

n 

calculated the sum of squares of residuals on p^, with p ,, taken to 
be p^. This latter sum was about twice as large as the corresponding 
sum using Eq. (67). We therefore conclude that the theoretical results 
are best represented by 
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2 2 

The values of the constants a^- a^ obtained from the least squares 
fit to Eq. (68) are given below: 

a^ = -0.0003089792906 
= 0.06065874406 
a^ = -0.000148501292 
a^ = 0.1740994651 

a^ = -0.01913423303 
a^ = - 0.0002122867568 

With Eq. (68), the rms deviation is 4.l4 x 10 The maximum abso- 
lute error on p calculated from Eq. (68) is 7.78 x 10 This error 
n 

is sufficiently small that the use of Eq. (68) is suitable for data 
analysis. When experimental calibration data are obtained, the 
question of using the forms in Eq. ( 63 ) or (66) and how P^Qp_^ should 
be calculated must be looked at again. 

Using the foregoing results we can now formulate the prediction 



-27- 


analysis problem for epithermal die-away in terms of specific quanti- 
ties. With Eq. (68) defining the relation between the variables 
y, and p^, the quantities defined in section I I -A are as follows: 


.1 


-1 

(69a) 

“V| 

(69b) 

-(Pj) i 

(69c) 

2 

(69d) 


(69e) 

-(Pj)? 

(69f) 

"=2"^Vi 

(69g) 

-aj -a^y, - 2a^(Pj), 

(69h) 


We will assume statistical weighting; i.e. w = 

1 

and w = — r— . Thus. 


2 ^ y 


w = 


1 

2 ^ 
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“ I 


(f ' a ) + Tf' '1 1 
W Y;' L Pp (Pp) J 
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(691) 
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k1 


= Z 




i = l 


/s 

L. 


(69j) 
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■ o' - 7 i ^ (C,)”- 

\ y Y' \ Pc P<r' U Li a, a, \ kl/ 


^ ^ k=l 1=1 ^ ' 


(69k) 


In evaluating Eq. ( 69 k); it should be emphasized that the quantities 

//^Vi ^ 

are calculated from the calibration data, the quantities F , 

^ /N /N 

F„ , F , s , and s refer to a measurement on an unknown, and we 
Pj ' y ’ Pj 

have assumed no correlation between the measurements of y and on 
the unknown (i.e. s^^ = 0). Also in Eq. ( 69 k), e is an estimate of 

all systematic errors which may be present. 

It would be desirable to make a study based on Eqs. (69a-k) and 
obtain predicted answers to the questions posed in the Introduction. 

An appropriate computer program, called PREDEPI, was written and suc- 
cessfully tested. For reasons which are based on the work described 
in the next three sections, no extensive calculations have been carried 
out with PREDEPI. However, one significant result emerged from some 



-29- 


preliminary calculations with the program. If either one or two 
different values of p are used in the calibration procedure^ i .e. 


only one or two different non-hydrogenous compositions^ the matrix C, , 
is i 1 1 -behaved and unphysical values of a are obtained. This indi- 
cates that the least squares data analysis would break down in such a 
case. Therefore, at least three different non-hydrogenous experimental 
compositions must be used for calibration. Information of this type 
which is available in advance can obviously save a considerable amount 
of time and expense when such information is used intelligently to 
plan an experiment. In fact, prediction analysis was developed for 
just this type of experiment planning. 


B. Experimental Studies 

1 . Tests on Location of He/Cd Detector in CPNE Probe 

One of the action items which resulted from the CPNE experi- 
menters' meeting at Goddard Space Flight Center on 19^ 20 November 19^9 
concerned flexibility in locating the ^He/Cd detector. It was recommended 
that some tests be carried out to answer this question. 

We have performed a series of epithermal die-away measurements 
on our basalt and crushed granite models. The experimental set-up is 
shown in Fig. 3. A Sandia neutron generator and tungsten attenuator 
were configured as shown. Although the location of a Nal detector is 
indicated in the figure, no actual detector was in place during these 

3 

tests. Four He/Cd detector positions were used as shown. Positions 
1 and 2 were on the model surface, and positions 3 and 4 were coaxial 
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wi th the generator. 

The data obtained on the crushed granite model are shown in 
Fig. 4. The effects we are looking for are easier to see for this 
model because it exhibits the slowest epithermal die-away rate of all 
our models. Data for positions 1, 2, and 3 are shown. The data for 
position 4 were very simiiar in behavior to those for position 3. The 
data for the three runs have been normalized at channel SO, or a delay 
of about 45 pseCj and the behavior of all three runs after channel 60 
is the same within statistics. 

From the results shown in Fig. k, we conclude that there is 
essentially no difference in the results which would be obtained in 
the various positions, provided that the die-away data at delay times 
greater than about 40 nsec are used in the data analysis. The effect 
of fast and high energy epithermal neutrons for the different positions 
is different and is clearly shown in the behavior of the peak in the 
time distribution (i.e., greatest for position 1, least for position 
3 or 4) . The neutron outputs per burst for the normalized data shown 
in Fig. 4 were in the ratio 1:1:3 for positions 1, 2, and 3, respectively. 

The question of the location of the He/Cd detector was thor- 
oughly discussed at an experimenters' meeting at Sandia Laboratories 
in Albuquerque on 21-24 June 1971. For reasons of economy of space 
it seemed desirable to place the detector somewhere within the symmetrical 
volume defined by the neutron generator, fast neutron attenuator, and Nal 
detector. After further study of possible probe configurations and 
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3 

available He detectors, it was recommended by C. A. McKeever of Sand i a 
Laboratories and J. A. Waggoner, now of Schlumberger Well Services, 
that a detector of 3 in. active length and in. active diameter 
be placed coaxially with the neutron generator tube, on the side 
opposite the neutron attenuator, and with its geometrical center ap- 
proximately 20.7 cm from the target. This location would presumably 
produce an epithermal die-away response intermediate between that in 
positions 2 and 3 in Fig. 4. The electronics for driving and controlling 
the neutron tube would be located in the region immediately above the 
•^He/Cd detector. 

2. Experimental Calibration 

( 2 ) 

As discussed in the Final Report on Phase II of our work, ' 
there remains the problem of satisfactory calibration of the epithermal 
die-away measurement as it is utilized in the CPNE. After considering 
several alternatives, we decided against having multiple water analyses 
carried out on samples of our rock models. There seemed to be little 
possibility that we would be able to adequately assess the errors involved. 

For calibration of the epithermal die-away measurement, it was 
decided to make a series of laboratory models which consist of stoichio- 
metrically known elemental compositions. Hydrogen content and the contents 
of all other elements would be determined from measured amounts of the 
constituents in the mixture. Hydrogen content can be easily varied by 
changing the amount of hydrogen-bearing material. In order to carry 
out this calibration program, several constraints must be met: 
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1) The models must be just large enough to be effectively 
semi-infinite in order to keep the costs of material 
down and to minimize the extensive amount of labor 
required to prepare the mixtures. 

2) The materials themselves should be available at a 
reasonable cost. 

3 ) The hydrogen content should span a range from essentially 
zero to the order of about 1%, and about ten values of 

Pu should be used. 

4 ) The hydrogen content must be almost all stoichiometric; 
absorption and adsorption of water vapor from the air 
should be smal 1 . 

5 ) The mixture of hydrogenous and non-hydrogenous materials 
must be homogeneous. 

6) As mentioned in section I I I -A, at least three different 
values of effective oxygen densities must be used for 
proper least squares analysis when Eq. (68) is used. 

Some experiments related to constraint 1 were carried out in which 

3 

the Sandia generator and He/Cd detector in position 2 (see Fig. 3 ) 
were used to obtain epithermal die-away data on the basalt and crushed 
granite models. Runs were taken with the generator and detector at 
various distances from the edge of each model. The distance at which 
the data become equivalent to those obtained in the center of the model 
was used as a guide in determining how large the hydrogenous model 
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mixtures must be. From an analysis of these data, it was concluded 
that the radius and depth of the hydrogenous models should be about 
50 cm or more. 

With constraints 2) -6) as a guide, we investigated several materials 
as to their suitability for hydrogenous and non-hydrogenous components 
in the proposed mixtures. Handbooks of material properties “and catalogs 
of prices were scanned; a material was a candidate if it was available 
at a reasonable price and in a strictly stoichiometric form, if its 
physical nature was such as to allow easy mixing with other materials, 
and if it did not absorb or adsorb water to any appreciable extent. 

After a considerable amount of literature searching, obtaining samples 
of possible material for testing, and carrying out some tests relative 
to the constraints, the following materials were selected; 

1. Pure silica ("frac" sand), 20-40 mesh 

2. Graphite, granular, 20-48 mesh. 

3. Sucrose ("sanding sugar"), ^22 ^11' mesh. 

Sanding sugar is coarser than ordinary cane sugar which is "extra fine" 
(it will pass through a 40-mesh sieve). Samples of each of the above 
materials were placed in shallow open containers and left for several 
weeks in the accelerator room where the die-away experiments would be 
carried out. The sample weights were determined before and after 
exposure to the ambient air, the relative humidity of which is controlled 

-4 

to be about 50%. The silica and graphite gained less than 4 x 10 % by 

weight; the sanding sugar gained approximately 0.06% by weight. Both 
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of these are acceptably low and can be readily determined for each 
model which is made. 

An experimental tank to contain the mixtures was ordered which 
was approximately 114 cm in diameter and 60 cm in height. The tank 
was mounted within a metal framework which was supported by a pair of 
jacks on casters. A track for the casters to run on was laid on the 
floor of the accelerator room so that the tank could be moved away 
from the accelerator target for filling and then moved back so that 
the target could be located 6.3 cm (~2.5 in.) above the surface of 
the filled tank and at any position along a diameter. 

In order to experimentally verify that (1) the Van de Graaff 
generator and Sandia neutron generator produce equivalent epithermal 
die-away data, and (2) the small tank model would act essentially like 
the large tank models, insofar as epithermal die-away is concerned, 
the small tank was filled with some spare Burnet Red crushed granite, 
the same kind of rock that is in our large outdoor crushed granite 
model. Several epithermal die-away runs were taken on both the small 
and large models. The source-detector positions referred to below 
are the ones shown in Fig. 3. Runs were obtained with both the Van 
de Graaff and the Sandia neutron generator on the small model; only 
the Sandia generator could be used on the large model. Statistics were 
moderately good, being 50-100 counts per 1-p,sec channel at a delay of 
50 p,sec. 

The following is a summary of the results obtained in the runs 
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described above: 

1) For detector position 4 on the small tank of granite, 
there is no evidence statistically of any difference 
between the Van de Graaff and Sandia generators. 

2) After a delay of about 30 |j,sec, there is no evidence 
statistically of any difference in detector positions 

2 and 4 with the Sandia generator. During the first 30 
psec the presence of the insulating oil in the generator 
is evident at position 2. 

3) Runs taken with the Sandia generator and detector in 
positions 2 and 4 on the large crushed granite model 
show almost the same behavior as the corresponding runs 
on the small model, although there is a si i ght hint of 
a statistical trend for the epithermal die-away rate to 
be faster on the latter. This could indicate that the 
experimental tests for minimum model size described at 
the beginning of section III-B-2 were not completely 
adequate in that the effects of vertical neutron leakage 
were not completely accounted for. 

Following the tests described above, the small tank was filled 
completely with the dry "frac” sand and a series of epithermal die- 
away runs was taken. A model consisting of only dry sand was chosen 
as the first one for experimentation because it was the easiest to 
prepare and because some problems which might exist with the selected 
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calibration procedure would materialize most prominently with a model 
containing no hydrogen. The experiments and results discussed below 
were obtained "stochastically”, in the sense that the nature and 
objective of each run were dictated by the results of preceding runs. 

A block diagram of the electronics and a schematic drawing of 
the experimental geometry are shown in Fig. 5. The cadmium-covered 

3 

He detector was placed on the surface of the material. The tritium 

target was located a distance above the surface which approximately 

corresponds to the position of the target in a Sandia generator when 

the latter is resting on the surface. The target-to-center-of-detector 

distance was arbitrarily chosen to be 15 cm. The configuration shown 

is easily reproducible on small and large models and with either the 

Van de Graaff or Sandia generator. The neutron detector electronics 

( 2 ) 

arrangement is straightforward and has been described previously . 

The master timing clock for the system was a crystal -control 1 ed Tektronix 
Model 180A Time Mark Generator. The clock pulses, which defined zero 
time for each cycle, were fed to a Rutherford Model All delay generator 
whose output was used to initiate time analysis in the TMC time analyzer. 
The purpose of delaying the start of time analysis is to permit measure- 
ment of the background immediately preceding each fast neutron burst. 
Clock pulses from the Time Mark Generator were also fed to an externally 
triggered variable frequency oscillator which, in this arrangement, 

served only as a pulse shaper to provide a pulse (the "initial pulse") 

2 + 3 

to switch the H beam onto the H target and to trigger a second delay 
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generator. The delayed pulse from this generator was used to switch 
the beam off the target. 

The results of the first run on the dry sand model are shown in 

Fig. 6. The Van de Graaff was pulsed at a rate of 2000 Hz with a 

burst width of 5 psec. A delay of 460 psec was used on the Model All; 

since the period between bursts was 500 psec, this allowed 40 psec of 

background to be analyzed. The time channel width was 1 psec. The 

data were accumulated for 4 x 10^ bursts over a 5.5 hour run. The 

average neutron output per burst was approximately 4520. The channels 

spanning the burst are plotted at the bottom as overflow; for these 

4 

the counts/channel scale should be multiplied by 10 . An average 
background of 395 counts/channel, calculated over the first 44 channels, 
has been subtracted from the count in each channel. 

The fact that the data were accumulated to such good statistics 
allows some features of the die-away curve, which normally are hidden 
in the statistical scatter, to be seen. The apparent behavior of the 
data in Fig. 6 is as follows. Early transients have disappeared by 
channel 80, which is about 30 psec after the end of the burst. From 
channels 80 to about l80 the die-away rate appears to be almost purely 
exponential, with a very slight hint of convexity toward the time axis. 
After channel 180 the rate of decay gradually decreases. This type of 
behavior is not expected on theoretical grounds. Calculated epithermal 
die-away curves given previously (c.f. Fig. 25 in Reference 2) exhibit 
a convex shape toward the time axis in the asymptotic late-time range. 
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In order to verify the expected theoretical behavior for dry 

sand specifically, a calculation of the response of an epithermal 

detector on a semi -1 nf 1 nl te medium of dry sand was carried out. This 

( 2 ) 

calculation utilized the theory developed previously and is the 
same type as described in section lll-A of the present report. This 
program, called THEOEPI, was run with the following dry sand parameters 


§ = 0.1087 

Y = 0.07661 

E = 0.1544 
0 = 2.243 
d = 4.601 

Values of r and 2 were 15 and 3.4 cm, respectively. The solid curve 
shown in Fig. 6 is the theoretically calculated response and has been 
normalized to the data in the region immediately following the burst. 
The slopes of the experimental and theoretical die-away curves are not 
too different out to about channel l40; however from there on the dis- 
agreement becomes catastrophic. The conclusion is that the data are 
probably being perturbed by some effect which is time dependent and 
which is distorting the shape. Distortion of this type usually comes 
from one or both of two possible sources; (1) tailing on the fast 
neutron burst, or (2) room return neutrons. 

On the next run made on the dry sand model, the Van de Graaff 
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generator was replaced with the Sandia generator. The same geometrical 
arrangement was used as is shown In Fig. S, and the same electronic 
system was used with the exception that the EP 110A and EP 210A were 
replaced by a Rutherford B7B pulse generator which was externally 
triggered by the Time Mark Generator. The results of this run are 
shown in Fig 7. The same experimental conditions were used as in the 
run shown in Fig. 6, except that the total number of bursts was 5 x 10^ 
and the neutron output per burst was about 285. When the data of Figs. 

6 and 7 are overlayed for comparison, it is obvious that they look 
quite different. The data taken with the Sandia generator show a much 
greater rate of epithermal neutron die-away than the corresponding run 
with the Van de Graaff. Also the data show a consistent trend of con- 
cavity toward the time axis. These differences should be considered 
in conjunction with a similar pair of runs, described previously, which 
were taken when the small tank was filled with Burnet Red crushed granite. 
In the crushed granite runs, there was no discernable difference between 
the two neutron sources. We attribute the large differences seen in 
Figs. 6 and 7 to the fact that the "frac" sand has essentially zero 
hydrogen content. The hydrogen present in the insulating oil in the 
Sandia generator plays a key role in determining the behavior of epi- 
thermal neutrons in the dry sand model. Although the crushed granite 
would seem to be relatively dry, we estimate that the total content of 
hydrogen in the small crushed granite model is the order of 10 times 
greater than that in the oil filling the Sandia generator. Thus, in 
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the small crushed granite model the hydrogen contribution from the 
insulating oil is of secondary importance, whereas in the dry sand model 
it is the dominant factor. This clearly indicates the importance of 
eliminating hydrogenous materials from the generator if die-away 
measurements on truly dry materials are to be made properly. 

In order to further elucidate the results obtained in the run 
shown in Fig. 6, a run was taken with the Van de Graaff target centered 

3 

over the sand tank and the He/Cd detector buried in the sand. The 
detector axis was vertical and displaced 2.5 cm parallel to the cylin- 
drical axis of symmetry of the model; the center of the detector was 
located 30 cm below the surface of the sand. The motivation behind a 
run of this type was to increase the ratio of the contribution to the 
counting rate from the model itself to that due to room-return neutrons. 

The results of the above run are shown in Fig. 8. Experimental 
conditions for this run were the same as for the data in Fig. 6, with 
the exception that a total of 3.S x 10^ neutron bursts were used at 
an average output per burst of approximately 3200. By comparison with 
Fig. 6, it is seen that the die-away rate is considerably faster with 
the detector buried than on the surface, there is a definite trend to 
convexity toward the time axis in the asymptotic region, and there is 
still evidence of a time-dependent background present. 

It was thought to be very doubtful that a large difference, such 
as that observed in Figs. 6 and 8, would exist in epithermal die-away 
rates between the surface and buried detector positions in a truly dry 
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medium. In an attempt to support this opinion, we made a theoretical 
calculation with program THEOEPl which was similar to the calculation 
described in connection with Fig. 6. The material parameters were the 
same as given before. Values of r and z were 2.5 and 35 cm, respectively. 
A value of z=35 cm corresponds to a physical depth of 30 cm because the 
former must include the extrapolation distance, d. The solid curve in 
Fig. 8 shows the theoretical result. Two points are significant. The 
disagreement between theory and experiment is much greater in Fig. 8 
than in Fig. 6, and the two theoretical curves in Figs. 6 and 8 
practically overlay. Thus, if the theoretical calculations are trust- 
worthy, at least in a relative sense, the experimental die-away rates 
should be little different for a detector located on the surface of 
the dry sand model and buried 30 cm into it. On the basis of the data 
and theoretical calculations in Figs. 6 and 8, we concluded that, (1) 
the surface detector run was badly distorted due to some effect or 
effects, (2) the buried detector run more nearly represents the true 
epithermal die-away in the dry sand model, although it too is being 
distorted somewhat, and (3) the THEOEPl program fails badly to account 
for the observed data. 

In order to investigate the problem of apparent asymptotic distortion 
of the epithermal die-away data, some diagnostic runs were taken. Figure 
9 shows the first of these, wherein the experimental conditions are the 
same as for Fig. 8 with the exceptions that the period was doubled to 
1000 p,sec, channel length was doubled to 2 psec, and a delay of 920 psec 
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was used. In this way, more of the time cycle could be observed. 

The data in Fig. 9 have not been corrected for background. There is 
clear evidence that, due to some effect, the behavior of the die-away 
at long delay times is reflecting the presence of a slowly decaying 
component. 

A run similar to that of Fig. 9 was carried out with the detector 
in the surface position. These data are shown in Fig. 10, The combin- 
ation of lower epithermal neutron flux at the model surface and reduced 
shielding of the detector by the model from room neutrons results in a 
lower ratio of signal to time-dependent background. 

As a final test for extraneous influences the sand was emptied 
from the tank, the detector placed in the "surface" position, and a run 
taken with a period of 2000 p,sec, channel width of 4 p.sec, and a delay 
of 1840 p,sec. The data are shown in Fig. 11. Two background components 
are evident. A time-independent background, with an average intensity 

of 311 counts/4 |isec channel, is present and is probably due to neutron 

• 2 + 
production throughout the source cycle by the H beam when it is de- 

3 

fleeted off of the H target. A time-dependent background is also present 
and probably results from neutrons produced at the target which travel 
to the walls of the accelerator room, make a few collisions, and return 
to the detector. The die-away of these neutrons in time will be strongly 
influenced by flight-time effects. For example, a 1-ev neutron requires 
about 250 psec to travel from the nearest point on the walls to the 
detector location. A flight time of 250 psec is 62.5 channels at 4 psec/ 
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channel . A pulse from such a neutron would occur at about channel 103 
in Fig. 11, and we see that this is essentially in the middle of the 
time-dependent decay. Thus the assumed source of the time-dependent 
background seems reasonable and is not inconsistent with the data. 

At this point two serious problems had presented themselves, viz., 
distortion of the data by a time-dependent background, which was presumed 
to be due to room return neutrons, and an apparent large discrepancy 
between experimental results and theoretical predictions. Since the 
proposed calibration scheme is based ultimately on theoretical arguments, 
it seemed necessary to solve the second problem before actively tackling 
the first. Also, the calibration experiments would be an irreversible 
process — once the materials were mixed and used, they could not be 
unmixed in order to repeat the experiments if something were found to 
be wrong with the basic approach. The next phase of this work involved 
an extensive theoretical study whose objective was to solve the second 
problem mentioned above. 

C . Theoretical Calculations Related to Experiments 

There are several possible sources of discrepancy between the THEOEPl- 
calculated die-away curves and the experimental data. It is possible that 
the basic theoretical model is not adequate. It is also possible that the 
finiteness of the dry sand model is playing a major role in determining 
the time behavior of the neutrons. This possibility is supported to 
some slight degree by the comparison runs on the large crushed granite 
model and on the small tank filled with crushed granite, as mentioned in 
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section I I 1 -B-2. 

As a first attempt to study the problem theoretically, some 
further calculations with program THEOEPI were carried out. Two 
different semi -inf ini te media were considered; these were dry sand 
and sand +0.2% hydrogen by weight. The material neutron parameters 
for dry sand were given previously. Those for sand +0.2 wt.% H are 
as follows: 

I = 0.2916 

Y = 0.7258 

S = 0.1940 
D = 2.060 

d = 3.662 

3 

The He/Cd detector response was calculated as a function of time for 
r = 2.5 cm and z = 3.4, 5, 10, 20, 30, 40, 60, 80, and 100 cm. The 
series of values of z was used in order to be able to plot vertical 
spatial profiles to see if the tank size might be influencing the measure- 
ments. Vertical profiles for a delay time of 50 psec are shown in Fig. 

12. The depth in a semi -inf ini te medium corresponding to the effective 
depth (including extrapolation distance) of the experimental tank is 
shown at z = 64.7 cm. It should be noted that the detector location 
used for the data given in Figs. 8 and 9 was 30 cm below the physical 
surface, or at a z-value of 34.7 cm in Fig. 12. 
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From the results shown in Fig. 12 it seems quite possible that the 
effects of finite model size could be quite large in the case of dry 
sand. The maximum of the vertical distribution in the semi-infinite 
case occurs at a z-value corresponding roughly to the bottom of the 
small tank. A boundary condition on the finite model is that the flux 
shall be zero at the extrapolated bottom of the tank. Thus, the shape 
of the vertical profile will be quite different for a semi-infinite 
medium consisting of dry sand and the small experimental tank filled 
with dry sand. This does not necessarily mean that there will be large 
differences in the epithermal die-away curves for these two cases, but 
such differences are strongly implied by the foregoing argument. Al- 
though no experimental results have been obtained for a small model 
consisting of sand +0.2 wt.% H, the vertical profile for this case 
shown in Fig. 12 would lead one to suspect that measurements made at 
depths shallower than about 30 cm might differ very little from those 
made on a semi-infinite medium of the same composition. 

Based on the results of the THEOEP I calculations shown in Fig. 12 
and the preceding discussion, it seemed necessary to theoretically in- 
vestigate more thoroughly the effects of model finiteness. It was 
deemed too difficult to try to develop a complete three-dimensional 
calculation at this time. As a reasonable compromise, which would 
nevertheless be informative, it was decided to develop two one-dimensional 
theories for the geometrical situations shown in Fig. 13. In one case 
a semi -inf ini te slab was taken to be in the region d < z < “, and the 
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extrapol ated plane where the flux is zero is at z=0. A plane source 
of fast neutrons is assumed to be at z=d. In the second case a finite 
slab occupies the region d < z < a — d. The extrapolated boundaries 
are at z=0 and z=a. A plane source of fast neutrons is assumed to be 
at z=d in this case also. In applying the results of the theoretical 
calculations, care must be taken to add the extrapolation lengths to 
the physical dimensions and to locate detector positions properly with 
respect to the boundaries in Fig. 13. 

Details of the solutions of these two problems are given in Ap- 
pendices 1 and 2. The work involved in obtaining the solutions was 
greatly facilitated by the mathematical similarities between the one- 
and three-dimensional problems. Appendices 1 and 2 are little more 
than additions to Appendix 3 in Reference 2, wherein details of the 
three-dimensional solution were given. The neutron fluxes for the one 
dimensional semi -i nf i n i te and finite slabs are given below: 


Semi - i nf in i te slab 


0(z,v, t) 


2 

TT 


^t 
- Y 


e 


(vDt)^ 



where 
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Finite slab 


0(z,v,t) 



gIVt 00 /nTT\2 

" I sr„(^) H„(v,t) (72) 

n=1 



The quantities appearing in the above expressions have the same 
meaning aS in Reference 2; their definitions are included in the 
List of Symboi s. 

Although the flux given by Eq. (70) or (72) is of interest itself, 
it is not the theoretical quantity which is to be compared with the 
results of die-away experiments. Instead we must calculate the theo- 
retical detector response R(z, t), where 

R(z, t) = J* d){z,v,t) e(v) dv (74) 

and e(v) is the speed-dependent He/Cd detector efficiency. It was 
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( 2 ) 

found previously that a useful and probably adequate analytic 
approximation to e(v) is 


e(v) 



, V < 1 cm/)isec 


cm/p^sec 


(75) 


Figure l4, which is reproduced from Reference 1 , shows the correct and 
approximate curves for e(v) . 

Two computer programs, called INFSLAB and F INSLAB, were written to 
evaluate the expressions given for the fluxes and to calculate the de- 
tector response by Eqs. (?4) and (75). The following characteristics 
of the programs are pertinent: 

(1) In INFSLAB, the upper limit on x in the integral in 
Eq. (70) was fixed at x=5 and the integral was eval- 
uated with a 1024 - point Simpson's rule. 

(2) In F INSLAB, the series was terminated when n reached 
a value such that 



(3) In both INFSLAB and FINSLAB, the integral for R(z, t) 



was evaluated by a trapezoidal approximation. 


( 4 ) Values of the Gamma function were obtained by routines 
developed by G. C. Wallick of this laboratory. These 
routines are based on rational polynomial approxima- 
tions and return values of F(x) which are correct to 
13 significant figures. 

A series of calculations was made with the two new programs and a 
modified version of the older THEOEP I program converted to a one- 
dimensional problem. The objectives of the various calculations varied^ 
since the answers to several questions were being sought. The calcula- 
tions and results are described below. 

The first series of calculations was with FINSLAB for dry sand 
and sand + 0.2 wt.% H. The neutron parameters for these two materials 
were calculated from Eqs. ( 42 ) - ( 45 ) and the constants given in Table 
1 . The bulk density of the sand was taken as 1.573 gm/cc (measured) 
and d was assumed to be , The material neutron parameters are 

Zj 

given below; 


Parameter 

Dry Sand 

Sand + 0.2 ' 

1 

0.1087 

0.2916 

S(cm"^ 

0.1544 

0.1940 

Y 

0.07661 

0.7258 

D(cm) 

2.243 

2.060 

d(cm) 

4.601 

3.662 

VQ(cm/p,sec) 

5200 

5200 
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Values of a=60 cm and z=5 and 30 cm were used. The following sets of 
values of t and v were used; 

t = 5,10,20,40,60^80, 100,120 |j,sec 

V = 1, 1.5,2,3,4,6,8, 10, 15,20,30,40,60,80, 100, 150,200 cm/p,sec 

For each combination possible among the variables, the calculations were 
carried out in two ways; (1) by using the convergence criterion given in 
(2) above, and (2) using 50 terms in the series for 0(z,v, t). The pur- 
pose was to verify that the proposed convergence criterion was satisfactory. 
It was found in every case that the two ways of calculating the flux 
agreed to at least 9 significant figures (the number printed out). The 
results also showed that the series expression for 0(z, v, t) in Eq. (72) 
is very rapidly convergent. In most cases only a few terms were required 
to meet the convergence criterion, and in the worst cases some 20 terms 
were needed. Rate of convergence of the series was considerably greater 
for dry sand than for sand +0.2 wt.% H. 

One result obtained from the above calculations was that the set of 
v-values had to be spaced finer in order to adequately outline a neutron 
velocity profile. In order to accomplish this, all further calculations 
with INFSLAB or FINSLAB were done with the following set of values of v; 
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Range of v 

Increment (Av) Betw 

1 < V < 2 

0.1 

2 < V < 5.2 

0.2 

5.2 < V < 6 

0.4 

6 < V < 10 

0.5 

10 < V < 20 

1 

20 < V < 50 

2 


The range of v's stops at 50 cm/p,sec. This upper limit is used^ not 
because the detector efficiency is negligible at this speed, but because 
the neutron fluxes are negligible for v > 50 cm/psec, even for the 
smallest value of t used (5 psec) . 

The next series of calculations was carried out for dry sand with 
both INFSLAB and FINSLAB. In each case the detector was assumed to be 
on the surface at z=d=4.6 cm. For FINSLAB the slab thickness was 60 cm, 
which gave a value of a=60+2d=69.2 cm. Epithermal die-away curves calcu- 
lated from the two programs are shown in Fig. 15 along with a repeat of 
the corresponding experimental run. The FINSLAB curve has been normalized 
to the data at the beginning of the die-away and the INFSLAB curve has 
been normalized to the FINSLAB curve in the asymptotic decay region. From 
the results in Fig. 15^ it is seen that the INFSLAB calculation and the 
data disagree badly over the entire time range, whereas the FINSLAB cal- 
culation and the data show rather good agreement over about 80 psec of 
the decay followed by an Increasing divergent trend. The agreement over 
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the early part of the decay is somewhat artificial, since it has already 
been pointed out that the experimental die-away is being held up by a 
time-dependent background. However, the significant point of Fig. 15 
is that the true decay can only be faster than is observed and thus 
the degree of agreement with the FINSLAB calculation and disagreement 
with the INFSLAB calculation must be greater than appears in the figure. 
This result reinforces the conclusion that the dry sand model is not 
effectively semi -inf ini te. 

Next, calculations were made with INFSLAB and FINSLAB for a model 
of sand +0.2 wt.% H. As in Fig. 15, the detector was assumed to be on 
the surface at z=d, the extrapolation distance, and the slab thickness 
was 60 cm for FINSLAB. Epithermal die-away curves for these two cases 
are shown in Fig. 16, together with the corresponding cases for dry 
sand shown in Fig. 15. The curves in Fig. 16 are normalized absolutely 
to each other, so that no vertical shifting was done in plotting. Several 
things are noteworthy. There is a persistent difference in both shape 
and intensity between the dry sand cases. The two cases of sand +0.2 wt,% 
H are very nearly the same, however, in both shape and intensity. This 
strongly suggests that a model of the size used in the experiments which 
consisted of sand +0.2 wt.% H would probably be a quite satisfactory 
model of the corresponding semi -i nf ini te medium. 

The next calculations that were carried out were for the purpose 
of comparison with the data in Fig. 8. The primary objective was to 
look into the possibility of a combination of experimental and theoretical 
calibration. If the theory, especially FINSLAB, could come reasonably 
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close to predicting the observed data, it would be a reasonable approach 
to calibrate a few theoretical parameters on the basis of some selected 
experiments. Calculations were carried out with both INFSLAB and 
F INSLAB for dry sand, a detector located at z=35 cm, and the standard 
set of values for t and v. For FINSLAB, the slab thickness was taken as 
60 cm (a = 69.2 cm). Figure 17 shows the calculated epithermal die- 
away curves and a repeat of the data given in Fig. 8. The disagreement 
between the data and the INFSLAB results is again quite large. The 
disagreement between the data and the FINSLAB results is not too bad 
when a rough estimate is made of the correction which should be removed 
due to the time-dependent background. This has not been done in the 
figure because there is no defensively objective way to carry it out. 
Even after a reasonable correction for the background, the negative 
logarithmic slope of the experimental data will be greater than that 
of the FINSLAB curve in the 40-60 p,sec delay region. This is not un- 
expected because the theory is only one-dimensional, whereas the experi- 
ment is three-dimensional and there is neutron leakage both radially 
and vertically. FINSLAB accounts for only the latter. The addition of 
radial leakage can only increase the epithermal neutron die-away rate, 
and it seems a reasonable supposition that a complete three-dimensional 
theoretical calculation might agree quite well with the observed die- 
away given in Fig. 17, provided that the time-dependent background is 
el imi nated . 

The next series of calculations which was carried out was related 
to the question of depth of Investigation. This problem has a direct 
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bearing on the sizes of models which would be necessary in order to 
look effectively semi -inf ini te for epithermal die-away measurements. 

The problem has an indirect bearing on the depths of investigation 
achieved in the thermal die-away and capture gamma-ray measurements, 
in the sense that the thermal neutron flux is driven by the epithermal 
flux. It should be stated at the outset, in regard to spatial questions, 
that the results of INFSLAB or FINSLAB should be considered as approx- 
imate guidelines rather than accurate predictions. This is because the 
assumption of constant neutron scattering cross sections, and the fact 
that cross section values in the epithermal range are used, is an ex- 
cellent one insofar as time dependence of the neutron flux is concerned, 
but a relatively poor one for spatial dependence. 

Calculations were carried out with INFSLAB for dry sand and sand 
+0,2 wt.% H, t=60 p,sec, the standard set of values of v, and thirteen 
values of z from 5 through 200 cm. Calculations were carried out with 
FINSLAB for dry sand and sand +0,2 wt.% H, t=60 p,sec, the standard set 
of values of v, and the following combinations of values of a and z: 

a = 70 cm; z = 5^10,20,30,40,60 cm 

a = 130 cm; z = 5, 10, 20, 30, 40, 60, 80, 100, 120 cm 

a = 210 cm; z = 5, 10, 20, 30, 40, 60, 80, 100, 1 20, 140, l60, 180, 200 cm 

The z- values used in INFSLAB are those above for a = 210 cm. A value 
of t = 60 p,sec was chosen as a representative delay time. The results 
of the calculations are shown in Figs. 18 and 19. In each case the 
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detector response, which is simply a certain velocity-weighted average 
of the flux, is plotted as a function of detector position, z. Note 

/f 

that the detector response has been multiplied by 10 for dry sand and 
10^ for sand +0.2 wt.% H. 

Several features of the results plotted in Figs. 18 and 19 are 
noteworthy. For dry sand and a=70 cm, the counting rates are very much 
less than those for a semi-infinite medium; in fact, the spatial dis- 
tribution is almost a fundamental mode, the antithesis of semi -i nf i ni te- 
ness. In order to have a proper model for epithermal die-away calibra- 
tion with the detector on the surface (z fa S cm), for dry sand the 
thickness should be greater than 70 cm but need not be greater than 
130 cm, whereas for sand +0.2 wt.% H a thickness of 70 cm is adequate. 

If epithermal calibration is desired for a buried detector position, 
primarily for corroborative and diagnostic purposes, larger models 
would be required, the sizes depending on the particular detector loca- 
tion. 

The results shown in Figs. I 8 and 19 are of interest in the impli- 
cations they contain for thermal die-away measurements by a gamma-ray 
detector and thermal neutron capture gamma-ray spectral measurements. 
Since the thermal neutron flux is derived from additional slowing down 
of the epithermal flux, and the detector responses in the figures were 
calculated at a delay of 60 p,sec (at this delay most of the response 
is due to neutrons with energy ~1ev), the spatial distribution of thermal 
neutrons and thermal neutron capture gamma rays must look something like 
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the curves given in Figs. l8 and 19. The most surprising implication 
is the degree to which a dry sand model of thickness about 70 cm fails 
to represent a sem i -i nf in i te medium. Considering the fact that the 
response of a Nal detector placed on the surface will be a volume 
weighted average over the neutron flux, the curves would indicate as 
much as an order of magnitude difference in counting rate between a 
truly semi-infinite dry medium and one whose thickness is of the order 
of 70 cm. 

Availability of all the INFSLAB and FINSLAB calculations which 

have been discussed thus far provided us with the opportunity to check 

the accuracy of the approximate expression obtained for the flux during 

( 2 ) 

the work of Phase II . The one-dimensional approximate solution is 
easily obtained by going through the three-dimensional solution and 
making the appropriate changes. The result for the flux is 



where 
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8 - vDt + H - *(^) - ij , (77) 

t(x) is the Digamma function, and all other symbols have the same 

meaning as given in the List of Symbols. The above flux is formally 

the same as that for the three-dimensional problem evaluated at r=0, 
with the sole exception that the exponent on 4np is 3/2 in the three- 
dimensional case instead of the 2 appearing in Eq. (76). 

The flux given by Eq. ( 76 ) was evaluated for dry sand and sand 
+0.2 wt.% H, for z=d, for t=5>60, and 120 psec for dry sand and 5^40, 
and 80 |J,sec for sand +0.2 wt.% H, and for the standard set of v 
values. Comparison with the exact fluxes in the corresponding 
INFSLAB calculations showed the following: (I) for dry sand the 

approximate flux is correct to within about 1%; (2) for sand +0.2 

wt.% H, the approximate flux is too small by 17% for t=5 p-sec, 12% 
for t=40 psec, and 8% for t=80 psec. Thus the slope of the approxi- 
mately calculated die-away curve should be almost exactly correct 
for dry sand, and have a negative logarithmic slope in the 40-80 psec 
range which is too small. By the law of propagation of errors, the 
fractional error in y, the slope defined by Eq. (36), is 



( 78 ) 
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assuming the curves to be normalized to be equal at t=40 p,sec and 
there are no errors in calculating R from 0 by Eq. (7^),. A calcula- 
tion for sand +0.2 wt.% H based on Eq. (78) gives |Ay/y 1 1.6%. 

Slopes calculated from the theoretical detector responses for dry 
sand and sand +0.2 wt.% H show discrepancies between THEOEPI and 
INFSLAB of 6.95 and 1.12%, respectively. The reason that the dry 
sand slope is off by about 7 % when it should be essentially correct, 
whereas the slope for sand +0.2 wt.% H is only off by about 1.1% 
when it should be off by 1.6%, almost surely lies in the fact that 
the detector response R is calculated in INFSLAB by a trapezoidal 
approximation. This is inferior to the 1024-point Simpson's rule 
used in THEOEPI. The discrepancy in the calculations is greater 
for dry sand than sand +0.2 wt.% H because of the shape of 0(z, v, t) 
as a function of v. Figure 20 shows plots of 0 versus v from INFSLAB 
for the values of the variables mentioned above. The fluxes change 
much slower with v for sand +0.2 wt.% H than for dry sand. Thus a 
trapezoid approximation for the area under the product of 0(z,v, t) 
and e(v) will be more accurate for the former case. 

D . Status of Analysis of Epithermal Die-Away Data 

Based on the results of the experimental and theoretical studies 
described in the two preceding sections, the status of the calibra- 
tion and analysis of the epithermal die-away measurement in the CPNE 
may be summarized in the following points: 
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1. A completely empirical calibration as originally conceived 
Is not possible with models which are k feet in diameter by 2 feet 
deep. Corrections for finiteness would have to be made in many of 
the cases^ and this would have to be based on theory. There still 
remains the problem of eliminating the time-dependent background if 
the models are used indoors. 

2. Construction of six (the minimum needed) or more models 
sufficiently large to be effectively semi -i nf i n i te is prohibitive, 
both from the cost and time standpoint. 

3. The best possibility for calibration in the future is 
probably a combination of experimental and theoretical approaches. 

The idea would be to use a few models outdoors where there would be no 
epithermal neutron return problem. The basic model might consist of 

a large tank of dry sand, roughly 400 cm in diameter by 200 cm deep. 

In the upper part of this tank, a smaller cylindrical volume would be 
maintained separate from the larger volume. The dimensions of this 
part might be about l40 cm diameter by 70 cm depth. When the small 
volume is filled with dry sand, a uniform dry sand model would be ob- 
tained which is effectively semi-infinite. Mixtures of sand and hy- 
drogenous material could be put in the small volume, the resulting 
models being also effectively semi-infinite. Cost of materials and 
time and labor of mixing would be minimized by this cyl inder-wi thin-a- 
cylinder arrangement. In these experiments it would be essential to 
use a generator with almost no hydrogenous material in it. This could 
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be achieved by use of the new Sandia generator if it is available^ or 
by repackaging an old one in a thin metal gas-tight housing and using 
an insulating gas, such as SF^. 

k. In order to utilize a greater portion of the epithermal die- 
away data for analysis, the possibility of getting a more convenient 
form of the theoretical solution should be considered. 

The question remains of what to use for an interim calibration 
so that the data obtained with the new test probe can be analyzed. 

The only possibility seems to be to use the theoretical approach given 
in section lll-A. This method is based on the use of Eqs, (68) and (6l) 
with the values of a^ through a^ given by the results of the least 
squares fit of Eq. (68) to the 27 theoretical cases. One way to get 
some idea of the adequacy of this approach is to analyze the data, 
shown in Fig. 2, taken on the four large rock models to see if the 
results seem to be at least reasonable. As a first step in carrying 
out this idea, the data were fit in the least squares sense over the 
time range 40 to 60 psec to the following formula: 

In R(t) = -yt + b (79) 

The values of y are the least squares estimates for the negative log- 
arithmic slope of the die-away curve. Values of for each model 
were calculated from the values of the nuclear constants given in 
Table 1, the elemental analyses given in Table 6 of Reference 2, and 
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the bulk densities of the models. Table 6 of Reference 2 is repro- 
duced in this report as Table 3- Finally values of were calculated 

H . 

from Eq. (68). These data are summarized in Table 4. The values of 
bulk density, p, are those obtained by accurate measurements as part 
of the thermal die-away studies to be discussed in section IV. The 
uncertainties on y are the standard deviation resulting from the fit 
of Eq. (79) to the data in each case. 

In order to reduce the results to a form which would allow com- 
parison with some independent measurement, the values of p„ were con- 

n 

verted to weight percent hydrogen and then to weight percent water. 
These were compared with the "best" chemical determinations of water 
content obtained previously. These best determinations were given 
in Table 12 of Reference 2 and are reproduced here as Table 5. Table 
6 gives a final summary of the results. One bit of editing has been 
done on the "best" chemical value of q, viz., that for Texas Rock- 
ville granite. The best value for total water content given in Table 
5 was judged to be the sum of the bound water determined by combustion 
at the uses and absorbed plus adsorbed water determined by weight loss 
from heating at our laboratory. But the solid block of Texas Rockville 
must have effectively almost zero absorbed and adsorbed water in the 
dense state when the epithermal die-away measurement was made. Most 
of the absorbed and adsorbed water reported must have accumulated on 
the sample after the rock was crushed for analysis. Therefore, the 
"best" chemical value of Wl, „ for Texas Rockville is taken to be the 
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bound water only. 

The agreement between the "best" chemical values of Wj^ q and 
those deduced from the die-away Is quite good and may be fortuitous 
to some extent. There Is no way for us to obtain values of un- 
certainties on w „ from the "best" chemical, primarily because we 
have no idea of the errors involved in the bound water analysis. 

The errors quoted on the values of w^ q from the die-away measurement 
are rough estimates only and have been obtained from Eqs. (28) and 
( 68 ): 


s ^ > (F s )^ + (F s )^ + 

Ph ~ ' y y' ' P5 


( 80 ) 


In writing Eq. (80), we have assumed that calibration errors are lumped 
into the systematic error e. From Eq. (68) we have 

^ “ -"2 "=5'=5 




For € we assume a value of 4. 1^x10 , the rms deviation quoted previously 

between correct values of p„ and those obtained from Eq. (68). Values 

H 

of a_-a^ were given in section 11 1 -A, and we arbitrarily assumed s to 
Zb 

be 5% of p_. The value of s was obtained from 


w 


H^O 
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As a point of interest it should be remarked that the reason we 

had so much difficulty previously in deciding how to proceed with the 

calibration process may be seen by comparing the values of y and w^ q 

{%, "best" chemical) in Table 6. It had been assumed that y and w . 

2 ° 

would be functionally related in a reasonably simple way. However 
it can be seen by comparing these values that this is not true. It 
was only after the theoretical arguments leading to Eq. (65) and its 
equivalent, Eq. ( 67 ), that a suitable relation between y and the 
material variables was found which ultimately permits a determination 
of w„ o from the die-away data. 
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IV. THERMAL DIE-AWAY 
A. Background Discussion 

Difficulties which have been encountered in analyzing thermal 

( 1 2 ) 

die-away data have been discussed previously. ^ ' These difficulties 
were of two types: (1) smoothing the data, and (2) extracting esti- 

mates of physically significant parameters from the data. A partial 
solution to the first type of problem was described in Reference 2 
and was implemented in a computer program called GUDFIT. Basically 
this program fits the logarithm of the observed counting rate to a 
polynomial in time in each of two or three time domains into which the 
thermal die-away data are partitioned. This procedure yields a good 
smooth approximation to the data in about 80% of the cases we tried. 

In the remaining 20% there was observed to be some type of undesirable 
feature in the smooth approximation. More to the point is the fact 
that the second type of problem still remained. This problem can be 
summarized as follows: according to available chemical analyses, there 

are essentially two models with values of macroscopic absorption cross 
section about 0.0065 cm \ and two with values about 0.012 cm how- 
ever, the observed thermal neutron and capture gamma-ray die-away rates 
indicated no such grouping. 

The foregoing difficulties led us to speculate that either the 

thermal die-away rates were much more sensitive to the thermal neutron 

diffusion coefficient, D, than we had previously assumed, or that the 

chemically deduced values of E were in error, or that both of these 
' a ’ 
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might be true. We also discussed the possibility of deducing the bulk 
density from the thermal die-away measurement if D were a sensitive 
parameter. 

Three items of future work were contemplated relative to the 
thermal die-away measurements in order to surmount these difficulties; 

3 

1. Obtain a new and independent set of both He and Nal die- 
away data using the Sandia neutron generator and the IITRI 
probe configuration. These would be run longer than usual 
in order to improve the counting statistics. 

2. Have additional chemical or elemental analyses carried out 
on samples of the models by one or more laboratories that 
could make some estimate of the accuracy of their results. 

3 . Attempt to empirically find a suitable functional form 
that would allow us to estimate values of physically sig- 
nificant parameters from the data. 

During the work of Phase III, all three of these have been essentially 
accompl i shed. 

B. Determination of T, Values For Rock Models 

As mentioned above, the calibration procedure for thermal die- 
away was complicated by the ambiguous results described in the Phase 
II report; namely, that four distinctly different lifetimes were ob- 
served for both neutron and gamma-ray measurements on the four rock 
models, whereas only two distinctly different values of macroscopic 
absorption cross section were computed from the compositions provided 
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by the wet chemical analyses. Thus If one assumes that the lifetime 
measurements are correct, then it follows that either unexpectedly 
large diffusional effects exist and contribute to the observed lifetimes, 
or the chemical analyses are incorrect. Several theoretical studies 
were carried out in an attempt to see if any reasonable theoretical 
model for the thermal die-away process could be found which would lead 
to strong diffusional effects. No such model could be found. Since 
no theoretical justification for the needed diffusional corrections 
could be demonstrated, it was supposed that either the chemical analyses 
erred in the amounts of major elements present in the rocks, or some 
strongly absorbing trace element(s) had been overlooked. A very good 
candidate in the latter category is the element gadolinium. It has a 
thermal neutron absorption cross section of 46,000 barns; hence only a 
few parts per million of this element are needed to cause significant 
changes in the value of for the rocks. The credibility of this 
explanation is reinforced by the fact that the average concentration of 
gadolinium in igneous rocks is reported to be approximately 7 parts 
per million . 

In view of the uncertainty as to the precise composition of our 

rock materials, it was concluded that the macroscopic absorption cross 

sections of our rock materials should be measured by an alternate 

pulsed neutron technique. Details of this method are given by Meadows 
/o\ 

and Whalen' ' ; only the principal features of the method will be discussed 
here. Basically, the procedure is to carry out thermal neutron die-away 
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measurements in a small container so that the fundamental spatial mode 
is dominant in the observed die-away rate. The method is absolutely 
calibrated by a measurement in pure water. The value of for an 
unknown is then determined from a measurement in a mixture of the un- 
known material and water. The observed lifetime^ t, is related to 
the thermal neutron diffusion and absorption parameters, D and E , of 

3 

the mixture by the equation 

^ = v(E^ + D B^) (81) 

where v is the most probable speed of thermal neutrons (2200 m/sec at 

the customary reference of 20°C), and B is an effective geometric 

buckling for the fundamental mode of the decaying thermal neutron 

population. The fundamental mode buckling is determined by the boundary 

conditions applied to the solution of the neutron diffusion equation. 

To what extent this parameter is "effective" will be discussed subse- 

2 

quently; for the present B is simply regarded as a known quantity from 
an initial experiment using only pure water. Since the diffusion param- 
eters of water are accurately known, it is possible to extract E^ of 
the rock from a lifetime measurement using Eq. (8l), provided D of the 
mixture is known. D can be computed with good accuracy from the 
chemical analyses for the major elements of the rocks and the known 
value of D for pure water. Trace elements are wholly unimportant 
because the water dominates the magnitude of D for the mixture. 
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A suitable reservoir is needed to contain the rock-water mixtures. 
Its size is important because large diffusional corrections are required 
if It is too small, and harmonic modes are persistent if It is too 
large. We used a 20x20x30-cm parallelepiped tank made of 2.5 nim thick 
Luci te sheets. It was covered completely with cadmium to screen out 
thermal neutrons returning to the reservoir from the room. The target 
of the Van de Graaff accelerator was positioned at the center of one 
20x20-cm face of the reservoir. Shielded access into the top of the 
reservoir was provided for a 6x6-mm cylindrical, enriched ^Lil detector 
and a 10-in. Lucite light pipe which optically coupled the detector to 
a photomultiplier tube. In all of the experiments reported here the 
detector was located at the center of the reservoir. Numerous measure- 
ments were made, however, with the detector located at the top of the 
rock-water mixtures in order to check the effect of locating the detector 
within the mixtures. The effect was always small, and agreed well with 
the magnitude predicted by diffusion theory. Centralization of the 
detector was preferred because much higher counting rates could be ob- 
tained, and the first harmonic modes are completely suppressed. Fast 
neutron bursts of 100 p.sec were used in every case in conjunction with 
a period between bursts of 3000 |j,sec. In all of the experiments dis- 
cussed in this section, the channel length of the TMC time analyzer 
was set at 16 [J.sec. 

The geometric buckling for the fundamental mode of the decaying 
thermal neutron population in a parallelepiped reservoir Is given by 
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= (-£^f * ^ fe)' ^ 

where a, b and c are the true physical dimensions of the reservoir, and 
d is the extrapolation distance which must be added to the physical 
dimensions to account for the fact that no thermal neutrons can return 
to the reservoir through the cadmium shielding. For the parallelepiped 
geometry d is related to the diffusion coefficient by the simple ex- 
pression, d=2.13D. When Eq. (82) is used to relate' the diffusion 
parameters with measured lifetimes for cases where the fundamental mode 
decay is slightly contaminated with higher modes, the buckling is re- 
ferred to as "effective". A measure of the contamination is provided 
by the extent to which the values of a, b and c must differ from their 

true values in order for Eq. (8l) to be satisfied. In the experiments 

2 

reported here this difference is very small. Note, however, that B 
varies from mixture to mixture through its dependence on the extrapola- 
tion distance d. 

Figure 21 shows the first die-away data taken for pure water. Ex- 
ponential decay with a single time constant appears to begin around 
channel 50. To be a little conservati ve, channels 55 and 85 were selected 
as the beginning and ending channels, respectively, for all lifetime 

computations. These computations were made using Behrens' method for 

(9) 

treating single-exponential decay . The lifetime of thermal neutrons 
In pure water was measured twice, the values obtained being 146.31 + *56 
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|j,sec, and 145.79 + .51 M’Sec. The average values of vD and vS given 

3 

by Martinho and Paiva^^^^ for four recent pul sed-neutron measurements 

in pure water (referred to 20°C) are 0.036137 + 0.000244 cm /psec and 

0.0048l8 + 0.000025 psec \ respectively. Using these values the 

effective buckling can be computed from Eq. (8l) to be 0. 056147 + 

2 

0.000938 cm . Since the height and width of the reservoir are identical, 

and equal to 2/3 of its length, we can let a=b=2c/3=l and solve for 1 

using Eq. (82). The result is 1=20.07 cm which agrees well with the 

true width and height of the tank. This value of 1 was used in all 

2 

subsequent computations of B for rock-water mixtures; corrections for 
variations in d through the diffusion coefficients were always made. 
Die-away measurements for the four different rock-water mixtures 
are shown in Figs. 22 through 25. Note that in every case the die- 
away is more rapid initially than for pure water, but decay with a 
single exponential time constant is established prior to channel 55. 

All of the crushed rocks used in these experiments, except the Texas 
Rockville granite, were used exactly as they were received from the 
supplier. No further crushing or sieving was needed. This was not 
true, of course, for the monolithic block of Texas Rockville granite. 

In our first attempt to measure its cross section, samples were chipped 
from the block, as randomly as possible, then crushed and sieved to 
obtain almost uniform chips of approximately 1/8" diam. ter. Temporarily 
set aside was approximately 2/3 of the original sample which had been 
reduced to a very fine particulate material by the crushing operation. 
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When die-away measurements were made on the water-saturated chips, and 
the data reduced to a value of S for the solid rock, it was discovered 
that the value of the cross section was too small to agree with die- 
away measurements made on the large block. This inconsistency was 
rectified, however, by including the fine particulate material in the 
lifetime measurement. To do this, it was necessary to saturate thin, 
alternating layers of rock chips and fine dust with distilled water. 

The fine material migrated significantly into the rock chips, but no 
gross separations occurred as frequently happens when one attempts 
to mix aggregates of two different sizes. Lifetime measurements were 
also made on two different mixtures of the fine material and water. 

These measurements confirmed that the fine part of the crushed rock 
did indeed have a larger cross section than the coarse component, the 
cross sections being 0 . 0226 ? cm ^ and 0.01823 cm respectively. 

One obvious explanation for this result is to suppose that a 
principal neutron absorber is concentrated near the grain boundaries 
of the primary minerals of the granite, perhaps in a mineral of small 
relative abundance. The crushing operation could have preferentially 
removed the secondary mineral leaving a "depleted" coarse component in 
the sieve. 

The absorption cross sections measured for the four igneous rocks 
are given in Table 7- Included also are the absorption cross sections 
calculated on the basis of the wet chemical analyses. Note that in 
every case the measured absorption cross section exceeds the cross section 
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computed using the chemical analyses. It would thus appear that at 
least one unknown thermal neutron absorber is present in all of the 
rocks, with the possible exception of the North Carolina Olivine 
(dunite). As mentioned previously, gadolinium is a prime suspect in 
the category of unrecognized offenders. The results of sufficiently 
sensitive specific chemical analyses for this element in our rock 
samples are not available. Two final points with regard to Table 7 
should be made. First, a small correction was applied to the measured 
cross sections to account for the presence of the detector and light 
pipe in the reservoir. The magnitude of this correction, in terms of 
the final cross section value, never exceeded 1%. Finally, the den- 
sities shown in Table 7 apply to the rocks as they are used in our 
large models. The density of the Texas Rockville granite therefore 
refers to solid rock, whereas the others are densities of crushed 
material . 

Solid-rock densities were needed for all of the rocks in order to 
calculate the absorption cross sections from the measured lifetimes. 
These were measured simply by adding 100 ml of crushed rock to a 200 ml 
flask containing 100 ml of pure ethyl alcohol and obtaining the weight 
of the added rock. Insolubility was assumed, but checked for each 
rock-alcohol mixture by measuring the density of the alcohol after each 
mixture was weighed. No measurable changes were observed. The results 
of these measurements are summarized in the following tabulation: 



Rock Description 


Density - gms/cc 


Balsam, N.C. Olivine 

3.24 

Burnet Red Granite 

2.64 

Knippa Basalt 

3.12 

Texas Rockville Granite 

2.63 


All of the figures cited for the densities in this tabulation are 
believed to be significant. 

C. New Thermal Neutron Die-Away Measurements on Large Models 

The configuration of the thermal neutron and gamma-ray probes 
used on the large rock models is shown in Fig. 26. These probes were 
always deployed on the surface of the models in such a way that the 
midpoint of the source-detector spacing coincided with the center of 
the model. In all of the experiments the repetition rate for the 
6-p,sec burst of fast neutrons from the Sandia generator was 200 pulses 
per second. A 25“Cm source-detector spacing was adopted for both the 
thermal neutron and the gamma-ray detectors. 

The first of the new die-away measurements on the large models 
was made with the neutron probe. An example of the data is shown in 
Fig. 27 as closed circles. In this example data were accumulated for 
2 hours on the crushed basalt model. Counting statistics are good, 
enough better than those obtained previously that the "tailing up" of 
the decay curve at large values of time is clearly evident. Such resul 
are frequently caused by thermal neutrons returning to an experiment 
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after scattering by a wall or ceiling. We had not anticipated problems 
of this nature because the models are located out-of-doors under a 
shelter of very light overhead structure. Nonetheless^ thermal neu- 
trons were being returned to the experiment by the 2"x6" ceiling 
joists and sheet iron roof, because a 1/32" sheet of cadmium completely 
covering the probe and the surface of the model eliminated the problem, 
as shown by the data in open circles in Fig. 27. It was subsequently 
verified that a single 6'x6' sheet of cadmium covering the experiment ' 
provided sufficient protection from the returning thermal neutrons. 

A complete set of die-away measurements using the cadmium-shielded 
neutron probe is shown in Fig. 28. The data are normal ized so as to 
display the differences in die-away rates observed for large values of 
time. Data accumulation times ranged from 1 to 3 hours, depending on 
the lifetime of the model. Adequate counting statistics were desired 
over at least 20 channels of the die-away spectrum where the decay is 
almost exponential. An attempt could then be made to fit this part 
of the data with the solution of the age-diffusion equation for the 
semi -inf ini te medium given in Reference 1; 


— ^ 4(0 + vDt)3 


n(r,z,t) = 


Se 


^4 tt( 0 + vDtjQ 




( 


-(.z-ll 


(Z-Hd)' 


4(0 + vDt) 4(0 + vDt) 

e — e 


(83) 


This equation was used in conjunction with the Gaushaus generalized 
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least-squares fitting routine' ' to f i t the data of interest shown 
in Fig. 28. In order to select, in some consistent way, the time at 
which the die-away period of interest begins, the data of Fig. 28 were 
first fitted by the empirical equation 


In 0(t) = 6^ “ ^2^ 


m 


The fitting began with channel 50 and stopped when less than 100 net 
counts were available. Very good fits were obtained, but the magni- 
tude of 6^ was uncertain typically by about 20%. It follows from Eq. 
(84) that 


O/at) In 0 = -e^ (85) 

Suppose we want the contribution of the second term in Eq. (85) to 
be a certain fraction of the first. Then 


263/1^02 = 3 , (86) 

where 3 'S the fraction. If channel 6l is judged to be the beginning 
of the asymptotic decay period for the Texas Rockville Granite, then 3 
can be calculated for this selection; the result is 3=.0?68. Using 
this value of 3^ the beginning channel for asymptotic decay was calcu- 
lated for the other materials with the following results. 
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Rock Description 

Channel Number 

Burnet Red Granite 

77 

Balsam, N.C. Olivine 

74 

Knippa Basal t 

64 


A1 though arbitrary, this procedure was believed to be better than an 
eyeball estimate for all four models. 

In the fit of Eq. ( 83 ) to the data shown in Fig, 28, the known 
values of D and Z were used, but the neutron age was treated as a 

3 

variable parameter. Table 8 gives the results of this study. Not 
unreasonable values were obtained for the ages, but the values were 
somewhat uncertain due to the statistical quality (and quantity) of 
the data. Attempts were made to fit larger portions of the data with 
Eq. ( 83 ); but poor results were obtained. This finding is not sur- 
prising because the simple theory cannot be expected to cope with the 
complicated slowing down and thermal ization processes which take place 
during the early portion of the die-away period. 

D. Method of Data Analysis 

However gratifying this limited success of age-diffusion theory 
might be from a theoretical point of view, the asymptotic fitting pro- 
cedure is basically unsatisfactory because a large amount of statis- 
tically reliable data is not used. An empirical remedy for this 
difficulty is Eq. (84). By comparison with Eq. ( 83 ) it can be seen 
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that 9^ should be proportional to the neutron source strength S, where- 
as 02 should be governed, at least for large t, by the magnitude of 
vS . All non-exponential behavior is clearly lumped into 9,. Hence, 

0^ can be viewed as a stripping parameter which allows more of the 
data to be used than would be possible with only a two-term expression. 

The die-away data shown in Fig. 28 were fitted with Eq. (84) over 
several ranges to check the sensitivity of 02 to the fit range. Best 
results were obtained when the time analysis commenced with channel 53. 


The same value of 02 was obtained as when either channel 52 or 54 was 
selected. These results are shown in the following tabulation of 02/v 
as a function of the beginning channel for the data fitting: 


Rock 

Ch 52 

Ch 53 

Ch 54 

Dunite 

0.00989 

0.00988+0.00010 

0.00986 

Burnet Red Granite 

0.01486 

0.01479+0.00023 

0.01474 

Basal t 

0.02006 

0.02003+0.00019 

0.02000 

Texas Rockville Granite 

0.02603 

0.02570+0.00066 

0.02619 


For beginning channels prior to channel 52 the value of was signifi- 

cantly influenced by the choice of the starting channel. The terminating 
channel of the die-away analysis period should be determined by statis- 
tical considerations. Data need not be rejected unless the scatter in- 
creases the uncertainty associated with the value determined for 02 . 

Examples of the data fits using Eq. (84) are shown in Figs. 29 through 
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32. Equation (84) gives slightly better fits to the entire set of data 
than 


In 0(t) = 0 '^ -Of^t + Q'^/t, 


(87) 


the second best of the simple forms. It was hoped that 0^ might be 
useful in determining the diffusion coefficient of the rocks, but this 
was not the case. For the data fits starting with channel 53 the statis- 
tical uncertainties associated with 0^ ranged from 20 to 50 percent. 

For a remote measurement, which could be allocated no more than one-half 
hour of data accumulation time, meaningful values of 0^ simply could 
not be obtained. 

It can be seen from Eqs. (83) and (84) that 0^/v should approach 
E in magnitude as E becomes very large. This trend is clearly evident 

d d 

in the following comparison of E^ and Q^/v: 


-1 


Rock 

Duni te 

Burnet Red Granite 
Basal t 

Texas Rockville Granite 


E — cm 
a 


0.00659+0.00030 

0.01149+0.00035 

0.01698+0.00030 

0.02113+0.00049 


G^/v — cm ^ 

0.00988+0.00010 

0.01479+0.00023 

0.02003+0.00019 

0 . 02570 + 0.00066 


A reasonable expression for correlating E and 0 /v would thus appear 

d 


to be 
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E /x = 1 — a,/x - a./x 
a 1 2 


a^/x' 


( 88 ) 


where 


X = e^/v 


(89) 


A calibration curve based on Eq, (88) is shown in Fig. 33. The co- 
efficients of Eq. (88) were determined by a least-squares fit to have 
the following values; 


. = 0.02555 

= - 0.001318 
a^ = 0.00002902 

An alternative form of Eq. (88) is clearly 

y = X — a^ — — a^/x^ 

where, for simplicity of notation, we have let 



(90) 


( 91 ) 
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According to Wolberg^^^, the variance of y is computed as follows. 
From the calibration data compute 


F ' = -Sy/Sa, = 1 

1 2 

F ' = -^y/Ba, = /x. 

3 I 

Form 


where 


+ y 


L. ~ O -riru f > 

I y. ' X X. ' 


= -By/9x. 


2 3 

— 1 — a,/x. — 2a«/x. 
1 1 ^ > 


(92) 

(93) 

(94) 


(95) 


(96) 


Calculate the matrix elements 


4 . . 

C, = E F 'F^'/L; 
kl 1=1 a, . 


(97) 


From a measured value of x on the unknown surface the variance of y is 
then given by 
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2 

y 


(F CT ) 

' X X 


z 


k=l 


k 

y F F c. , 

Lt a. a, kl ’ 

1=1 ' 


(98) 


2 

where a is the variance of the measured value of x, 

X 

The entire computational problem of fitting the calibration data 
using Eq, (88)^ and calculating S from a measured value of 0^/v is 
carried out by a computer program named SIGMA. The value of 
its fractional standard deviation are obtained from GAUSHAUS when the 
die-away data are fitted using Eq. (84). When the calibration data 
were processed by program SIGMA as unknowns^ the predicted cross sec- 
tions compared favorably with the measured values as they should: 


Rock 


Duni te 

Burnet Red Granite 
Basal t 

Texas Rockville Granite 


Measured T, 

0.00659+0.00030 

0.01149+0.00035 

0.01698+0.00030 

0.02113+0.00049 


Predicted E 

0 

0.00658+0.00035 

0.01167+0.00043 

0.01642+0.00033 

0.02160+0.00083 


The predicted fractional standard deviations of E vary from case to 

a 

case but average 3-7%. While greater precision is certainly desirable, 
such results are believed to be adequate and consistent with the com- 
bined experiment. Note that the calibration range is 0.010 < < 0.026 

in units of cm \ Since absorption can be expected to dominate the 
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measured lifetime for values of G^/v exceeding 0.026, extrapolation on 
the high side is probably safe enough. Diffusional effects are clearly 
important for O^/v < 0.01, hence extrapolation on the low side is 
probably very risky. 

It should be noted finally with regard to thermal neutron exper- 
iments that die-away measurements were made with the helium-3 detector 
buried 60 cm below the surface of the rock models. The detector was 
located on the axis of the cylindrical model directly below the Sandia 
generator which was resting on the model surface. These measurements 
were made to check the age-diffusion theory prediction that the life- 
time should increase only slightly with depth, and to make sure that 
the return neutron problem had been solved by the cadmium shielding. 
Figure 3^ illustrates a typical result of these experiments. The data 
are not normalized; except for small variations in the output of the 
neutron generator the two-hour measurements are directly comparable. 

Very similar asymptotic die-away rates are plainly evident; lifetime 
values calculated for the two cases differ by only 6 percent. In 
order to make this measurement on the block of Texas Rockville Granite, 
it was necessary to drill a I- 5 ” diameter horizontal access hole 60 cm 
below the surface of the slab. Similar results to those shown in Fig. 

34 were obtained. 

E. New Capture Gamma-Ray Die-Away Measurements 

Capture gamma-ray die-away measurements on the large rock models 
essentially duplicated the experiments performed with the thermal neutron 
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detector. No changes In timing or accelerator operation were made. 

The gamma-ray detector system was always biased to reject photons 
having energies less than 250 kev. Although it was not necessary to 
shield the gamma-ray detector from return neutrons, suitable shielding 
for the model surface would have been desirable. Unfortunately, com- 
pletely satisfactory shielding was either too expensive, or judged to 
be sufficiently thick that its use would interfere with the experiments. 
As a result none was used. it vould appear from the data that the 
effect of return neutrons is small if not negligible, quite in contrast 
to the situation encountered when thermal neutrons were detected. This 
is due to the fact that the gamma-ray detector responds to thermal 
neutron captures occurring at considerable depth within the model. 

Figure 35 is the gamma-ray counterpart of Fig. 28. Four distinctly 
different die-away rates are plainly evident, in agreement with the 
results obtained using the helium-3 detector. Fits of the individual 
die-away curves using Eq. (84) are shown in Figs. 3& through 39. As 
before, Eq. (84) gave the best fits to the entire set of data. The 
data fitting always began with channel 44. Earlier beginning channels 
could not be used because the rapid die-away component present in 
channels before 43 would clearly contribute to the analysis. Values 
of 02 /v obtained when the time analysis commenced with channels 43, 44, 
and 45 are shown in the following tabulation; 
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Rock 

Ch 43 

Ch 44 

Ch 45 

Dun i te 

0.00976 

0.00975+0.00020 

0.00970 

Burnet Red Granite 

0.01348 

0.01355+0.00056 

0.01318 

Basal t 

0.02008 

0.02007+0.00021 

0.01985 

Texas Rockville Granite 

0.02518 

0.02351+0.00102 

0.02280 


The time-base stability of the fitting procedure beginning with channel 
44 is clearly quite good. Note, however, that the statistical errors 
are somewhat larger here than for the neutron detector. 

The gamma ray calibration curve based on Eq. ( 88 ) is shown in Fig. 
40 . Least-squares values for the coefficients are 

a^ = 0.02452 

= -0.001699 

a^ = 0.00004246 

Using the computer program referred to earlier for calculating the magni- 
tude and variance of E given 0 -/v and its fractional standard deviation, 

3 ^ 

the following results were obtained for the gamma-ray calibration data 
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Rock 


Duni te 

Burnet Red Granite 
Basal t 

Texas Rockville Granite 


Measured 
a 

0.00659+0.00030 

0.01149+0.00035 

0.01698+0.00030 

0.02113+0.00049 


Predicted £ 
a 

0.00660+0.00056 

0.01138+0.00085 

0.01752+0.00037 

0.02063+0.00108 


The fractional standard deviations calculated by program SIGMA average 

5.8% for the gamma-ray data. Although significantly larger than the 

3.7% average obtained for the thermal neutron detector, these results 

do not appear to warrant giving serious consideration to the inclusion 

of a thermal neutron detector in the combined experiment. As with the 

thermal neutron calibration, extrapolation should be safe enough for 

G^/v exceeding 0.025 cm but not for values less than 0.01. 

A remote surface with less than 0.01 cm ^ could be a rock 

composed largely of quartz and containing no strongly absorbing trace 

elements. This case could probably be handled by adding a calibration 

point for clean, quartz sand. We note in this connection that the 

macroscopic absorption cross section of a clean quartz sand was measured 

recently by the technique outlined earlier. The value obtained for S } 

a 

0.00280 cm \ agreed almost exactly with the cross section computed on 
the basis of the chemical analysis provided by the supplier. 

A more serious problem is the possibility of encountering a very 
low-density rock having a very small macroscopic absorption cross section. 
Such a material could have a measured value of 02/v falling within the 
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calibration range, but one whose die-away rate would be governed by 
neutron leakage rather than absorption. 
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V. ANALYSIS OF VENUSIAN ATMOSPHERE AND 
SURFACE BY LIGHTWEIGHT CPNE PROBE 

A CPNE experimenters' meeting was held at I IT Research Institute 
in Chicago on 7 October 1970 to discuss and organize a presentation to 
be made to planetology personnel in OSSA. The primary purpose of the 
meeting and presentation was to consider and propose the use of a very 
small lightweight CPNE probe on a surface mission to Venus. One new 
idea brought out was the possibility of utilizing the CPNE to analyze 
the dense Venusian atmosphere during the rather slow descent of the 
spacecraft through the atmosphere. The presentation was made at NASA 
headquarters on 22 October. As a result of the presentation, and dis- 
cussions during and following the presentation, it was decided that 
high priority be given to studying the capabilities of a CPNE probe 
weighing less than 10 lbs. 

A special problem which must be looked into is that of analyzing 
the Venusian atmosphere. The basic question is which, if any, of the 
measurements in the CPNE could yield useful information when carried 
out singly or repeatedly during the descent through the atmosphere. In 
order to answer this question one must estimate counting rates due to 
the desired signals and due to sources of interfering background. As 
was pointed out in the Introduction, this ultimately involves predictions 
of the variances of the variables being determined. 

Data analysis procedures for the epithermal and thermal die-away 
measurements made on semi -inf ini te igneous rock models have been discussed 
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in preceding sections. From these results it should be possible to 
make some predictions about the capabilities of these two measurements 
in the Venusian atmosphere. Preliminary elemental detection thres- 
holds have been obtained for semi-inf ini te models for the inelastic, 

(12 3 13 ) 

capture, and cyclic activation gamma-ray spectral measurements. ’ ’ 

No similar estimates have been made for the case of an infinite medium 
such as the atmosphere of Venus, the composition of which is primarily 
CO^ plus small amounts of N2 and 

We were asked to investigate the capabilities of the CPNE for 
Venusian atmosphere analysis on an "as-time-perm i ts" basis. If the 
CPNE does prove to have significant promise for atmospheric analysis, 
it would considerably strengthen its position as a candidate for the 
first soft-landing mission on Venus. Unfortunately, time and money 
limitations did not permit us to look very extensively into this problem. 
This section describes the results of those studies which we were able 
to carry out. 

It is unlikely that any approach could result in absolute predictions 
of detection thresholds or parameter variances which could be relied on 
as to accuracy. For this reason, it was decided to attempt to obtain 
relative values between the infinite Venusian atmosphere and semi-infin- 
ite rock models for which we have data. If this could be done, the 
performance of the CPNE in the atmosphere of Venus under given experi- 
mental conditions could be related to its performance in a known situation. 

We consider first the epithermal die-away measurement. Recall that 
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the detector response is proportional to R(r,z, t), where 


R(r,z, t) 


J* 0(r,z,v, t) e(v) dv , 


( 99 ) 


0(r,ZjV,t) is the velocity- and time-dependent epithermal neutron flux 

at a detector located at (r,z), and e(v) is the relative detection 

efficiency. An approximate expression for the flux for the case of a 

( 2 ) 

point source on a semi -i nf ini te medium was given previously. If 
we assume that time channels of length t and a total of N bursts are 
usedj the number of counts accumulated in time channel i will be pro- 
portional to Y., where 


Y. 

I 



R(r,z, t) dt 


Time channel i 


psTINtR ( r, z, t . ) , 

where T] is the effective geometrical efficiency of the neutron detector. 
Using the expressions previously given for the flux and e(v), the number 


of counts in channel i for a semi -inf i ni te medium becomes 
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(1-0.715 — ) X 
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rn 


dx 


§Et. 

Y 


( 100 ) 


where 




Define the integral in Eq. (100) to Jj,^> so that 


1.43 HtNO 

3/o 2 

(4TT)'2r(i) 


j! 

I 


( 101 ) 


In Eq. (100) we have assumed that the detector is on the physical surface 
(z=d) . 

The expression for the number of counts in channel i for an infinite 
medium, such as the Venusian atmosphere, is similar to Eq. (100) except 
that the infinite medium flux must be used. This flux may be obtained 
from an analysis similar to that used in Appendix 3 of Reference 2. The 
approximate expression for the flux is 
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0(r,2^v, t) 



r Y ” 43 


V2 


(4tt3) 


( 102 ) 


In analogy with Eq. (101) we have, for the counts accumulated in channel 
i for an infinite medium. 


where 


Y I 1.43 TItNO j 


( 103 ) 


(4tt) ^(^) 


J.' = 

I 


00 

I („. 




' -X -’^2 ■ '•6 

e dx , 

(104) 


715 e 3 


yx 


and we have assumed that the source-detector distance is r, the same 
as for the semi -inf i ni te medium. 

Although Eqs. (101) and (IO 3 ) may not be extremely accurate in 
an absolute sense, the ratio of counts for the two situations is prob- 
ably given to reasonable accuracy by 
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Htnq 

n-) 
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( 105 ) 


In order to estimate the variance of a hydrogen determination 
from an epithermal die-away measurement in the Venusian atmosphere^ 
we use the prediction analysis version of Eq. (80): 




/S/A. 
2 2 
F a 

y y 


2 2 


A 

2 


+ F + e 


( 106 ) 


2 2 2 

Predicted values of F , F , and e must be obtained from assumptions 

y Pj' 

about the nature of the atmosphere and from a determination of the 


em . 


constants a,— a^ in Eq. (68) appropriate to the infinite medium probl 
I b 

2 

A predicted value of a must come from an analysis similar to the 

present one for the gamma-ray spectral portions of the CPNE. The method 

2 

of predicting values of a is discussed below. 

y 

A general treatment of the prediction of variance for situations 

for which an exponential decay model is valid is given in Chap. 6 of 

(4) 

the book by Wolberg. We have assumed that the epithermal die-away 
data will be analyzed over the range 40 to 60 |j,sec by using the form 



Y(t) 


Ae 


-yt 


(107) 
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and least squares fitting 1 nY to a linear relation in t to determine 
y. For the above form, Wolberg shows that the predicted values of the 
matrix elements, defined in Eq. (10), are as follows: 



I 






^12 “ ^21 



At.e 

I 


r2yti 




n 


I 


1,2 2 -2yt. 

A t.e ' I 



( 108 ) 


( 109 ) 


( 110 ) 


In the above equations it has been assumed that there is no error in 

the independent variable, t., and that it is measured in channels rather 

than time. Since 1-p,sec channels are usually used for the epithermal 

die-away measurement, the numerical value of y wi 1 1 be the same, whether 

measured in (p.sec) ^ or (channels) The quantity n is the number of 

2 

data points used and Oy is the predicted variance on the accumulated 

i 

count in channel i representing the die-away signal. 

Since there may be a non-negl i gi bl e time-independent background 
present (e.g. from an RTG), we assume that the observed count in channel 
i is C., where 

I 
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C. = Y. + b, 

I I 


( 111 ) 


and b i s an average background per channel determined from the portion 
of the data immediately preceding the fast neutron burst. This quan- 
tity would be determined by averaging the observed counts in some 
number, say m, of channels; 


m 



j=l 


( 112 ) 


Now we have 



and by the law of propagation of errors, 



But 


( 113 ) 



and 
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C. 

J 
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Therefore, 
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pioa^ 






(1 
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Equation (Il4) is to be used in Eqs. (108) - (110). 

XN /A, 

The sums in and Z^^ are difficuit to evaluate. How- 

ever Wolberg points out that, if the independent variable values are 
evenly spaced, to a sufficiently good approximation we can write 


C 


n 


= I 


A 

-2yt. 


1=1 Ae^^' + (1 + ^) b 


(115) 


_2yt 


A 1 /s 

+ (1 + -) b, 


Ae 


average 


( 116 ) 


w 


b a 


-2yt 


Ae"^^ + (1 + — ) b 
' m' 


dt , (117) 


where t and t. are given by 
a D 


*^a ~ ^1 2 


( 118 ) 


^b = ^ ^ 2 


(119) 
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In Eqs. (118) and (119)> t is the channel length. If the convention 
Is made that, insofar as the least squares fitting of Eq. (107) to 
the data is concerned, the first channel used in the fit will be channel 
1, and we assume that l-jisec channels will be used, we have 


A 

C 


n+0.5 

r 


11 


A 

-2yt 


V/ 

0.5 


A J^t 1 A 

+(! +-) b 
' m 


dt 


( 120 ) 


Similarly, the following approximate expressions 
be obtained; 


/A 

for ^22 


^2 ^ 


n+0.5 

-A 

V 

0.5 


te 


-2yt 


A , A 

Ae^ + (1 + —) b 
m 


dt 


( 121 ) 


n+0.5 


'22 




0 




A 

-2yt 


dt 


A 

Ae'T" . 


m 


( 122 ) 


A predicted estimate of o can be obtained from the matrix elements 
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^ir ^12^ ^22 follows: 



/\ 

c 

2 

^ -1 

11 

a *: 

y 

^22 /N 

/s 


^11 * ^22 ^12 


( 123 ) 


In principle, an estimate of the variance of a hydrogen determin- 
ation from an epithermal die-away measurement in the Venusian atmos- 
phere can be made by using Eqs. (105), (106), (120), (121), (122), and 
(123). Although time did not permit us to carry out any numerical 
calculations, we feel that the foregoing analysis is reasonable and 
could provide meaningful results if and when it seems desirable and 
justifiable to obtain them. 

Analyses similar to the foregoing could be carried out for the 
thermal die-away and gamma-ray spectral measurements. These analyses 
will probably be more complicated than the epithermal die-away problem, 
but at the least the above approach might serve as a model upon which 
the other analyses could be based. Active work in this area will have 
to await the decisions required to commit the necessary funds and man- 
power to the task. 


s,-. 
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VI . MISCELLANEOUS TOPICS 
A. Neutron Output Monitor 

The possibility and desirability of remotely measuring the neutron 
output of the neutron generator which will be used in the CPNE has been 
discussed informally in the past among the experimenters. Such a 
measurement could be useful in several ways. It would provide a positive 
indication of how the generator was performing. A comparison between 
observed counting rates in the parts of the CPNE and the measured 
neutron output could yield some additional information about the surface 
being analyzed that is not available otherwise. Finally, normalization 
of measurements made in the CPNE to a known neutron output would allow 
direct comparisons among several sites which were analyzed In a mission 
that employed some type of roving vehicle or similar means for carrying 
out multi-site analyses. 

In laboratory measurements with the Sandia generator, the neutron 

output has been measured by either activation of Cu foils through the 

Texas Convention^ or activation of F in Teflon cylinders, followed 

( 2 ) 

by counting in a Nal well crystal. Clearly such methods are impossible 
for remote applications of the CPNE. A satisfactory monitoring method 
must be an integral part of the CPNE probe and must be capable of oper- 
ating within the overall constraints of the experiment. 

We have done some preliminary tests on a monitoring method which 

3 

would utilize the He/Cd detector used for the epithermal die-away measure- 
ment. The only additional items required to make the monitoring measure- 
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ment would be a pulse height discriminator and a scaler. The method is 

based on detecting 1^-Mev neutrons produced by the generator by means 
3 3 3 

of He recoils from the He(n,n) He elastic scattering reaction. From 
the mechanics of elastic scattering, the maximum energy which would be 
imparted by a neutron with energy E to a nucleus of mass M (relative 
to the neutron mass) is 

E«,ax = [' - (^)'J ^ 

For M=3 and E=l4 Mev, this gives E =10.5 Mev. The reaction by which 

^ ^ max ' 

3 3 

epithermal neutrons are detected is He(n,p)'^H, which has a Q-value of 
0.77 Since epithermal neutrons bring very little kinetic 

3 

energy to the reaction, the pulses produced in the He detector from 
them are much smaller than those from the ^He recoils due to l4-Mev 
neutrons. Pulses from the recoils will extend from 10.5 Mev down to 
zero, and some reasonable bias level, corresponding to a lower energy 
value, must be selected in order to achieve adequate counting rates. 

The question that then must be answered is whether this lower bias 
level will permit fast neutrons which have been scattered off the sur- 
face being analyzed (and possibly off of various structural parts of 
the probe) to be detected in sufficient quantity as to cause the pro- 
posed monitoring scheme to be too environment sensitive. 

In order to investigate the possibilities of realizing the method 
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3 

proposed above, the He/Cd detector used in the epithermal die-away 
measurements described in section IH-B was set up near the Van de 
Graaff target in a position identical to that used for the epithermal 
experiments. Initially the small tank was empty. A pulse height 

. 3 

spectrum of pulses from the He/Cd detector was taken with the Van de 
Graaff producing l4-Mev neutrons in a steady state condition. This 
spectrum is shown as solid points in Fig. 4l . For comparison purposes, 
an Am-Be source was put in a paraffin moderator and the spectrum of 
pulses recorded. These data are shown as open circles in Fig. 4l . The 
spectra have been normalized to equal intensity at channel 20. 

The spectrum from the moderated Am-Be source shows the familiar 
epithermal distribution rising to a peak at channel 20, which corres- 
ponds to an energy of 0.77 Mev, and then dropping rapidly to a compon- 
ent which tails off to higher energies and presumably results from 
detection of fast neutrons which have escaped from the moderator. The 
spectrum from the D-T neutron source shows an intense continuum which 
monotonical ly decreases towards higher energy. It was decided to use a 
bias level corresponding to 3 Mev to carry out the ensuing monitor 
tests. 

3 

With the He/Cd detector positioned over the empty tank, a series 
of runs was made at various beam currents on the Van de Graaff. For 
each run, a Teflon cylinder in our standard monitoring position was 
exposed and counted in the Nal well crystal to measure the neutron 
output. Following this, the small tank was filled with crushed granite 
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3 

and the He/Cd detector repositioned exactly the same as in the runs 
with the empty tank. A series of runs was made at various beam 
currents and the neutron output monitored by Teflon activation. 

All observed count rates were corrected for background and each 
of the two sets of data were fitted to the form 

3 1 3 

■^He recoils = k( F counts). 

The values of k were as follows: 

Empty tank: k = 0.96 + 0.04 

Granite filled: k « 1.15 + 0*06 

There thus seems to be a statistically significant environmental effect 

3 

on the high energy He-recoil monitor. It is in the expected direction, 

since the presence of the granite will scatter some fast neutrons into 

the He/Cd detector and increase the He recoil counts, thus increasing 

the observed value of k. Since the effect is only about 20%, this might 

be an acceptable influence. The high-energy bias point could also be 

3 

raised to diminish the environmental effect, but at the expense of He 
recoil counting rate. 

The standard conditions which were used in these tests were as 
follows: 1" active diameter by 4" active length, 10 atm. ^He/Cd detector; 

cylindrical axis of detector passed through '^H target at 0° to the H 
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beam; geometrical center of detector located 11.4 cm. from the target. 

Under these conditions^ we observed a monitoring sensitivity of 6.75 
5 3 

xIO neutrons/ He recoil count. For the anticipated combinations of 

repetition rates and neutron output per burst, we estimate that the 

3 

above sensitivity would lead to approximately 100 He recoil counts 
per minute. Unless the ambient fast neutron background is unusually 
high, this should be an acceptable counting rate. The possibility of 
using a monitoring technique based on the foregoing approach should be 
considered for the future. 

• • 

3 

B. Flight-Qualified He Detector 

One problem which we had not specifically looked into during the 

3 

Phase I and II work is the availability of flight-qualified He detec- 
tors of the type needed to do the epithermal die-away measurement. As 
a result of an inquiry made in July 1970, It was learned that such 
detectors are available from Reuter-Stokes of Cleveland, Ohio. A 
series of Reuter-Stokes' Model RS-P4-0806-10 detectors were tested by 
Reuter-Stokes under contract to TRW, Inc., and passed all requirements 
for space flight. This detector model has a 1" active diameter and an 

3 

active length of up to 8". He fill pressures up to 20 atm. are 
aval 1 able. 

As a result of the CPNE meeting at Sandia Laboratories in June, 

3 

1971^ it appears that a 1" diam. x 6 or 8" active length He detector 
would not fit into the new CPNE probe which is being designed as well 
as a detector which has a shorter active length and a larger diameter. 



3 

By keeping the volume of He gas approximately the same, however, it 
should be possible to keep the epithermal neutron counting rate un- 
changed . 

A second inquiry was made of Reuter-Stokes, and they stated that 

a Model SK -678 could be modified to meet our requirements. Two versions 

were suggested: a I 2 in. diam. x 3 in. length counter or a in. diam. 

X 4^ in. length counter. Either one could be filled with 10 atm. of 
3 

He plus a small amount of and would operate at less than 2000 volts. 

It appears that the shorter detector is probably most suitable from the 
standpoint of dimensional considerations in the new probe. 
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VII. SUMMARY 

It Is evident from the nature of this report that our main occupa- 
tion during the work of Phase III has been in the area of data analysis 
and interpretation of the epithermal and thermal neutron die-away measure- 
ments. By and large, our efforts have been successful. Although it 
would have been desirable to carry out an experimental calibration of 
the epithermal die-away measurement, we feel that the theoretically- 
based data analysis procedure which has been devised is adequate for 
the immediate application at hand, viz., the analysis of data which will 
be obtained with the new lightweight CPNE probe. The successful analysis 
of the thermal neutron and capture gamma-ray die-away data, and the un- 
ravelling of the problems associated with it, is particularly gratifying 
since this situation had represented something of a "thorn-in-the-side" 
matter to us since the work done in Phase I. Analysis of thermal neu- 
tron capture gamma-ray die-away data obtained with the new CPNE probe 
in the upcoming field test will be carried out by the technique described 
in the present report. 

Some problems associated with the die-away measurements remain 
which need working on, but these are not regarded as having sufficiently 
high priority to justify proposing another phase of work in the near 
future. Rather, these can be reconsidered at the time, if and when, the 
CPNE is taken off the shelf to get it prepared for a campaign to have 
it accepted for use on a suitable approved mission. Specific problems 
which could be worked on are as follows: 
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1. The basic method of analyzing the epithermal neutron die- 
away data should be extended, probably on the basis of theo- 
retical work, in order to permit more of the data obtained 
in a measurement to be used in the analysis. 

2. The finite theory of epithermal die-away must be worked 
out for a true three-dimensional cylinder before a combined 
experimental and theoretical calibration can be satisfac- 
torily carried out. 

3. Actual calculations involving prediction of variance should 
be carried out for both epithermal and thermal die-away for 
cases of future interest. 
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LIST OF SYMBOLS 


Y 


i 




observed value of the dependent variable at the i ' th data point 
standard deviation of Y. . 

I 

observed value of the j ' th independent variable at the i ' th 
data point. 



standard deviation of X. . 

J • 

calculated (or adjusted) value of the dependent variable at 
the i ' th data point. 



calculated (or adjusted) value of the j ' th independent variable 
at the i ' th data point. 





true (or mean) value of the dependent variable at the 1 ' th 
data point. 

true (or mean) value of the j'th independent variable at the 
i 'th data point. 


n 

m 



P 


= number of data points. 

= number of independent variables, 

= calculated values of the unknown parameters, 
= true values of the unknown parameters. 

= standard deviation of a. . 

k 

= unbiased estimate of a 
= number of unknown parameters. 
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neutron speed. 

speed of source neutrons. 

neutron lethargy. 

macroscopic scattering cross section 
diffusion coefficient. 

average logarithmic energy decrement per collision. 

characteristic constant defined in equation (^1). 

microscopic scattering cross section. 

mass of nucleus/mass of neutron.. 

weight fraction of element. 

extrapolation distance. 

gamma function, 
digamma function. 
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APPENDIX 1 


We wish to solve the Boltzmann transport equation for the semi- 

infinite slab shown in the upper part of Fig. 13. This problem is 

( 2 ) 

similar to the three-dimensional problem solved previously and the 
solution outlined below will lean heavily on the exposition given in 
Appendix 3 of Reference 2. 

The same assumptions and boundary conditions are used in the 
one-dimensional problem as in the three-dimensional one. With deflni 
tions and notation the same as used before, the equation to be solved 
along with the appropriate boundary conditions is as follows: 


^e^|f+E0-D^ = 

Ot s ^ 

0 oz 


sj ® ^ 6(u'-u)J 0 du' + 6(z— d) 6(u) 6(t) 

0(0, u, t) = 0(z, u, 0) = 0 
lim 0(z,u, t) = ' ® 

Z “ 2 -t 00 

In Eq. (A1.1) we have explicitly assumed a unit plane impulse source 1 
cated at z=d. The second part of the Integral in Eq. (Al.l) can be 


(A1.1) 


A1-1 



carried out, and the equation to be solved becomes 




at 


Sz 


JL5 

2 

y 


! 


u '~u 

u Y 

e 0 du' + 6(z— d) 6(u) 6(t) 


(A1.2) 


Define the Laplace transform of 0 with respect to t and the 
Fourier sine transform with respect to z by 


0^^(k,u,s) = J sin kzdz J 


® 0(*‘^ z, u, t) dt 


Taking the Laplace - Fourier transform of Eq. (A1.2), and using the 
boundary and initial conditions and well-known properties of the trans- 
forms, we obtain 




u u 

^ 0*''' du' 


+ 



sin kd 6(u) 


Define a new function 


iLF . 

^ by 


(A1.3) 



(A1.4) 


) 


Al-2 



Substituting 0*"^ from Eq. (A1.4) 


into Eq. (A1 .3) gives 


iJ 

(is! ^ ^ „,2V UF ^ 

Vvq Y P 



(A1.5) 


The next series of steps in the solution is formally identical 
with those following Eq. (A3. 6) in Appendix 3 of Reference 2. Noting 
that 'll and 0 differ only at zero lethargy and recalling the fact that 

3 

the He/Cd detector response is essentially negligible at u=0, we 
have the analog of Eq. (A3. 9) of Reference 2: 

2 - 

0(z,v,t) = n e ^ l(z,v,t) (A1.6) 


where 


l(z,v, t) = 


00 

dk sin kz sin kd 

0 


2 

-vDtk 

e 


g(k) 


and 
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g(i<) = 



In obtaining Eq. (A1.6) we have changed variables from lethargy to 
velocity by using the relation 


«5(z,v, t) = 0 (z,u, t) 1^1 

It is convenient to change the variable of integration in the expres- 

2 2 

sion for l(z, v, t) to x = vDtk . When this is done we have finally 


_ 

2 e 

t) = - T 

"" (vDt)2 



where 
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It should be remarked, as a point of comparison, that the closed 

form expression for 0(r,z,v, t) which is given in Eq. (A3. 12) in Refer- 

— f (x) 

ence 2 was obtained by writing g(x) as e ' '', expanding f(x) in a 

2 

Taylor's series about x=0, and retaining terms through x . The accu- 
racy of this approximation was unknown and, consequently, the useful- 
ness of the closed form solution remained to be demonstrated. We 
avoid the above approximation in this treatment and leave 0(z,v, t) in 
the form given by Eq. (A1 .7) . 
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APPENDIX 2 


In this appendix we wish to solve the Boltzmann transport 
equation for the finite slab shown in the lower part of Fig. 13. 
The basic equation to be solved is Eq. (A1.1). The boundary and 
initial conditions are 

0(0, u, t) = 0(a,u, t) = 0 

0(z, u,0) = 0. 

Define the Laplace transform of 0 with respect to t by 


0^(Z,U,S) = J 


e 0(z, u, t) dt 


Taking the Laplace transform of Eq. 
condition we obtain 



(A1.2) and using the initial 


Lf ^ 

J e ^ 0*" du ' + 6(u) 6(z— d) 

0 


In order to satisfy the boundary conditions at z=0 and z=a, we 
try a Fourier series solution of the following form: 


(A2.1) 


A2-1 



0'"(z,u, s) = ^ 0^(u,s) sin 


n=l 


It is also necessary to expand 6(z— d) in a Fourier series: 


= I A„ (^) 


(A2.2) 


n=1 


By a straightforward calculation, we find 


2 . / nnd \ 

A = - sin I ) 

n a \ a / 


Substitution of the above quantities in Eq. (A2.1) yields the 
following equation which must be satisfied by 0|^(u, s): 


+ "(t) = 


u'— u 


J e 0n(uSs)du' + | 


Y 0 


(A2.3) 


Define a new function \{f (u, s) by 

n 


A2-2 



2 


(A2.4) 


s) = 0^(u, s) 




2 

- sinl . 
a \ a / 


5(u) 


s_ 

V, 




“( t ; 


Substituting 0 (u^s) from Eq. (A2.4) into Eq. (A2.3) gives the following 
n 

equation to be solved for ([f (u, s): 

n 




ilf^(u, s) = 


IS 

2 

Y 


u '-u 



ilf^(u', s) 


u 



(A2.5) 


From this point on in the development, as in Appendix 1, the series 
of steps is formally the same as those following Eq. (A3. 6) in Appendix 
3 of Reference 2. Again noting that 'jr^ and 0^ differ only at zero 

3 

lethargy where the He/Cd detector response is very small, the analog 
of Eq. (A1.6) in the preceding appendix is 


0(z,v, t) = - e 

a 


_ lIYt «> 
n=l 


i nTT \ 

1 


— vDt \ 3 J 


h (v, t) (A2.6) 
n 


A2-3 



where 






2/y 

Dv/nTT\^ 

+ §E\. aJ 


(A2.7) 


Similar to the final result given by Eqs. (A1.7) and (A1.8) in 
Appendix 1, no approximations have been made in arriving at Eqs. (A2.6) 
and (A2.7) that are not reasonable and of an estimable accuracy. 


A2-4 



TABLE 1 - VALUES OF NEUTRON-NUCLEUS CONSTANTS 


Element 

M. 

1 

4 


CT. (Barns) 

r . 
1 

r. 

t 

r 

0 

r. 1 - 
1 ^1 

0^0 

r . V. 
1 1 

r 5 Y 
0^0 '0 

H 

1 

1 .000 

1.000 

21 

12.64 

- 

- 

- 

C 

12 

0. 158 

0. 108 

4.7 

0.236 

1 .65 

2.17 

2.87 

0 

16 

0. 120 

0.0816 

3.8 

0. 143 

1.00 

1.00 

1.00 

Na 

23 

0.0845 

0.0571 

3.3 

0.0864 

0.604 

0.426 

0.298 

Mg 

24 

0.0811 

0.0548 

3.5 

0.0878 

0.614 

0.415 

0.279 

A 1 

27 

0.0723 

0.0488 

1.4 

0.0312 

0.218 

0.131 

0.0783 

Si 

28 

0.0698 

0.0470 

2.2 

0.0473 

0.331 

0.193 

0.111 

K 

39 

0.0504 

0.0339 

2.0 

0.0309 

0.216 

0.0907 

0.0377 

Ca 

40 

0.0492 

0.0330 

3.0 

0.0452 

0.316 

0. 130 

0.0524 

Fe 

56 

0.0353 

0.0237 

11.5 

0.124 

0.867 

0.255 

0.0741 




Rock 

Type 


Basal t 


Grani te 


Dun i te 


TABLE 2 - SUMMARY OF THEORETICAL EPITHERMAL 
DIE-AWAY CALCULATIONS 


p 


Ph 



% 

y 

1.5 

0 

0 

1.0425 

0.84180 

0.75420 

0.01329 

2 

0 

0 

1.3900 

1.1224 

1 .0056 

0.01663 

2.5 

0 

0 

1.7375 

1.4030 

1.2570 

0.02379 

1.5 

0.0015 

0.00225 

1.0425 

0.84180 

0.75420 

0.05288 

2 

0.0015 

0.00300 

1 .3900 

1.1224 

1 .0056 

0.07274 

2.5 

0.0015 

0.00375 

1.7375 

1.4030 

1.2570 

0.09310 

1.5 

0.0030 

0.00450 

1.0425 

0.84180 

0.75420 

0.08473 

2 

0.0030 

0.00600 

1.3900 

1.1224 

1.0056 

0.11446 

2.5 

0.0030 

0.00750 

1.7375 

1 .4030 

1.2570 

0. 14420 

1.5 

0 

0 

1.0251 

0.90645 

0.82650 

0.01390 

2 

0 

0 

1 .3668 

1.2086 

1 . 1020 

0.01820 

2.5 

0 

0 

1 . 7085 

1.5108 

1.3775 

0.02883 

1.5 

0.0015 

0.00225 

1.0251 

0.90645 

0.82650 

0.05491 

2 

0.0015 

0.00300 

1 .3668 

1.2086 

1 . 1020 

0.07581 

2.5 

0.0015 

0.00375 

1 . 7085 

1.5108 

1.3775 

0.09723 

1.5 

0.0030 

0.00450 

1.0251 

0.90645 

0.82650 

0.08687 

2 

0.0030 

0.00600 

1.3668 

1 .2086 

1.1020 

0.11752 

2.5 

0.0030 

0.00750 

1.7085 

1.5108 

1.3775 

0. 14820 

1.5 

0 

0 

1.1224 

0.93360 

0.83325 

0.01408 

2 

0 

0 

1 .4966 

1.2448 

1.1110 

0.01885 

2.5 

0 

0 

1 .870B 

1.5560 

1 .3888 

A A-S 1 1 -7 
V • 1 1 / 

1.5 

0.0015 

0.00225 

1.1224 

0.93360 

0.83325 

0.05572 

2 

0.0015 

0.00300 

1 .4966 

1.2448 

1.1110 

0.07699 

2.5 

0.0015 

0.00375 

1.8708 

1.5560 

1 .3888 

0.09882 

1.5 

0.0030 

0.00450 

1 . 1224 

0.93360 

0.83325 

0.08768 

2 

0,0030 

0.00600 

1 .4966 

1 .2448 

1.1110 

0. 1 1861 

2.5 

0.0030 

0.00750 

1 .8708 

1.5560 

1 .3888 

0.14959 


All densities in gm/cc 


y in units of |j,sec 




TABLE 3 - ELEMENTAL COMPOSITIONS OF ROCK MODELS 
(WEIGHT PERCENT) 



Burnet Red Granite 
P = 1.38 gm/cm^ 

Texas Rockville Granite 
P =2.63 gm/cm^ 

Knippa Basalt 
P = 1.60 gm/cm 

Dunite, Balsam, N 
P = 1.75 gm/cm^ 

0 

47.74 - 1.06 

47.82 - .76 

43.95 - 2.51 

44.12 - .97 

Si 

33.06 - .75 

32.66 - .45 

19.13 - 1.69 

20.37 - .27 

A 1 

7^79 

8.11 

5.53 

0.52 

Fe 

2.24 

2.78 

9.07 

7.04 

Ca 

0.94 

1 .30 

9.09 

0.24 

Mg 

0. 14 

0.06 

8.53 

27.19 

Na 

3.02 

2.95 

2.09 

0.01 

K 

3.57 

3.77 

1 .20 

0.008 

Ti 

0.20 

0.25 

1.92 

0.0006 

Li 

0.001 

0.001 

0.0017 

< 0.00005 

B 

< 0.0003 

< 0.0003 

< 0.0003 

- 

Cl 

0.01 

0.01 

0.01 

0.02 

1 

1 

P 

0,04 

0.07 

0.37 

0.01 

H 

0.077 

0.075 

0.316 

0.077 



TABLE h - SUMMARY OF ANALYSIS OF EPITHERMAL 


DIE-AWAY DATA 


Rock Type 

,p (gm/cc) 

p^(gm/cc) 

y(p,sec ’) 

p^(gm/cc) 

Burnet Red 
Grani te 

1 .38 

0.7935 

0.0323 + 0.0037 

0.001090 

Texas Rockville 
Gran i te 

2.63 

1.517 

0.0606 + 0.0029 

0.001533 

Knippa Basalt 

1 .60 

0.9024 

0.0942 + 0.0067 

0.005017 

Dun i te 

1.75 

1.072 

0.0485 + 0.0041 

0.001645 




TABLE 5 - WEIGHT PERCENT WATER IN ROCK MODELS 



USGS 

12-6-67 

Mobil FRL 
(NMR) 

Mobil FRL 
(weight loss) 

USGS (Bound) 

Mobil FRL (absorbed) 

Burnet Red Granite 

Absorbed H 2 O 

0.08 

- 

0.13 

0.13 

Bound H 2 O 

0.56 

- 


0.56 

Total H 2 O 

0.64 

0.36 

0.49 

0.69 

Texas Rockville Granite 

Absorbed H 2 O 

0.05 

- 

0.12 

0.12 

Bound H 2 O 

0.55 

- 

0.19 

0.55 

Total H 2 O 

0.60 

0.37 

0.31 

0.67 

Knippa Basalt 

Absorbed H 2 O 

0.57 

- 

0.82 

0.82 

Bound H 2 O 

2.0 

- 

1 .42 

2.0 

Total H 2 O 

2.57 

0.92 

2.24 

2.82 

Dunite, Balsam, N.C. 
Absorbed H 2 O 

0.05 


0.03 

0.03 

Bound H 2 O 

1.0 

- 

0.61 

1 .0 

Total H 2 O 

1 .05 

- 

0.64 

1.03 





Rock Type 

Burnet Red Grani 
Texas Rockville 
Knippa Basalt 



TABLE 6 - COMPARISON OF ANALYSES OF WATER CONTENT OF ROCK 
MODELS DEDUCED FROM DIE-AWAY AND FROM COMBUSTION 



y(|isec"^) 

rt( 7 o,"BEST" CHEMICAL) 

H2O 

w {%, DIE-AWAY) 

H2 0 

te 

0.0323 + 0.0037 

0.69 

0.71 ± 0.14 

Grani te 

0.0606 + 0.0029 

0.55 

0.52 + 0.07 


0.09^2 + 0.0067 

2.82 

2.83 + 0.28 


0.0485 + 0.0041 

1 .03 

0.85 + 0.13 


Balsam Dunite 



TABLE 7 - MACROSCOPIC ABSORPTION CROSS SECTIONS FOR ROCK MODELS 


Model 

^a(CHEM. ANAL,)-cm“’ 

^a (MEASURED)*-cm"^ 

Dens i ty-qm/cc 

Balsam, N.C. Olivine 

0.00592 

0.00659 + 0.00030 

1.75 

Burnet Red Granite 

0.00601 

0.01 149 + 0.00035 

1.38 

Knippa Basalt 

0.01 160 

0.01698 + 0.00030 

1 .60 


Texas Rockville Granite 


0.01215 


0.02113 + 0.00049 


2.63 



TABLE 8 - NEUTRON AGES DETERMINED BY 
LEAST-SQUARES FITS TO THERMAL DIE-AWAY DATA 


Mode 1 

D-cm 

V -1 

*^a - cm 

0 - cm^ 

Fit Range - 

Channels 

Balsam, N.C., Olivine 

1.55 

0.00659 

676 + 30 

74 

- 140 

Burnet Red Granite 

2.17 

0.01 1A9 

1244 + 179 

77 

- 120 

Knippa Basalt 

1.33 

0.01698 

638 + 42 

64 

- 108 


Texas Rockville Granite 


1.13 


0.02113 


297 + 40 


61 - 90 



(Y., X.) OBSERVED POINT 

I I 



\ 


\ 
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FIGURE I 

TYPICAL DATA ANALYSIS SITUATION (TWO-DIMENSIONAL) 

SHOWING RELATION BETWEEN OBSERVED, CALCULATED, AND 

TRUE DATA POINTS. R AND R ARE RESIDUALS; E AND E ARE ERRORS. 
yi Xi ^yi Xi 





COUNTS / 1/iSEC / 30 MIN 


FIGURE 2 

EXPERIMENTAL EPITHERMAL DIE-AWAY DATA 


He /Cd DETECTOR 

SANDIA GENERATOR IN IITRI PROBE 
2500 PULSES / SEC 
OUTPUT ~ 400 NEUTRONS / BURST 
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FIGURE 3 

BLOCK DIAGRA/Vl OF ELECTRONICS FOR 
EPITHERMAL DIE-AWAY MEASUREMENTS 













FIGURE 6 





THEOEPI CALCULATION 
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FIGURE 8 



THEOEPI CALCULATION 
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FIGURE 9 

EPITHERMAL DIE-AWAY ON DRY SAND MODEL 
DETECTOR BURIED 30 CM 

VAN DE GRAAFF 
lOOO-nSEC PERIOD 
2 iiSEC/CHANNEL 
2.5 X 10^ BURSTS 
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FIGURE 10 

EPITHERMAL DIE-AWAY ON DRY SAND MODEL 
DETECTOR ON SURFACE 

VAN DE GRAAFF 

iooo-m,sec period 

2 p,SEC/CHANNEL 
2.5 X 10^ BURSTS 




FIGURE 11 


EPITHERMAL DIE-AWAY WITH EMPTY TANK 
DETECTOR IN SURFACE POSITION 
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VAN DE GRAAFF 
2000-nSEC PERIOD 
4 nSEC/CHANNEL 
1 . 1 X 10^ BURSTS 



RELATIVE DETECTOR RESPONSE 


FIGURE 12 
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FIGURE 13 

GEOMETRICAL CONFIGURATIONS FOR 
ONE-DIMENSIONAL CALCULATIONS 


He /Cd DETECTOR EFFICIENCY 









He/Cd. DETECTOR RESPONSE X 10 


FIGURE 16 






FIGURE 17 
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VERTICAL PROFILES OF ^He/Cd DETECTOR 
RESPONSE FOR SAND + 0.2 WT. % H 

t = 60 fj.SEC 
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FIGURE 24 

THERMAL DIE-AWAY OF CRUSHED BASALT 
SATURATED WITH PURE WATER IN SMALL 
RESERVOIR 
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CHANNEL NUMBER 






FIGURE 26 


PROBE CONFIGURATION FOR 
THERMAL DIE-AWAY EXPERIMENTS 



NEUTRON PATH 













FIGURE 28 

THERMAL NEUTRON DIE-AWAY DATA FOR 
FOR ALL FOUR LARGE ROCK MODELS 
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FIGURE 29 


THERMAL NEUTRON DIE-AWAY DATA FO R 
DUNITE MODEL FITTED O VER CHANNELS 
50 - 140 


= 9 ,- 6 z 1 + 63/ 1 
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FIGURE 34 

COMPARISON OF THERMAL DIE-AWAY ON 



CHANNEL NUMBER 
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FIGURE 36 

GAMA/IA-RAY DIE-AWAY DATA FOR DUNITE 
MODEL FinED OVER CHANNELS 44 - 90 

In (f> (t) = ^1 - ^ 2 ^ + ^ 3 /t^ 
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